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CONVERSION  FACTORS,  U.S.  CUSTOMARY  TO  METRIC  (SI) 
UNITS  OF  MEASUREMENT 


U.S.  customary  units  of  measurement  used  in  this  report  can  be  converted 
to  metric  (SI)  units  as  follows: 


Multiply 

by 

To  obtain 

inches 

25. 4 

millimeters 

2.54 

centimeters 

square  inches 

6.452 

square  centimeters 

cubic  inches 

16.39 

cubic  centimeters 

feet 

30.48 

centimeters 

0.3048 

meters 

square  feet 

0.0929 

square  meters 

cubic  feet 

0.0283 

cubic  meters 

yards 

0.9144 

meters 

square  yards 

0.836 

square  meters 

cubic  yards 

0.7646 

cubic  meters 

miles 

1.6093 

kilometers 

square  miles 

259.0 

hectares 

knots 

1.852 

kilometers  per  hour 

acres 

0.4047 

hectares 

foot-pounds 

1.3558 

newton  meters 

millibars 

1.0197  x  10-3 

kilograms  per  square  centimeter 

ounces 

28.35 

grams 

pounds 

453.6 

grams 

0.4536 

kilograms 

ton,  long 

1.0160 

metric  tons 

ton,  short 

0.9072 

metric  tons 

degrees  (angle) 

0.01745 

radians 

Fahrenheit  degrees 

5/9 

Celsius  degrees  or  Kelvins1 

JTo  obtain  Celsius  (C)  temperature  readings  from  Fahrenheit  (F)  readings, 
use  formula:  C  -  (S/9)  (F  -32). 

To  obtain  Kelvin  (K)  readings,  use  formula:  K  =  (5/9)^  (F  -32)  ♦  273.15. 


8 


SYMBOLS  AND  DEFINITIONS 

d  water  depth 

f  wave  frequency 

fHF  high-frequency  spectral  cutoff  used  for  confutations 

fLF  low-frequency  spectral  cutoff  used  for  confutations 

fp  frequency  at  the  middle  of  the  spectral  band  of  maximum  energy 
density 

g  acceleration  due  to  gravity 

H  wave  height 

Hs  significant  wave  height 

L  wavelength 

Lp  wavelength  corresponding  to  a  wave  of  frequency  f p 
m^  nth  moment  of  the  spectrum 

P(h^)  probability  associated  with 

Qp  spectral-peakedness  parameter 

<\n  nth  moment  of  the  sea-surface  elevation  distribution  function 

S(f)  spectral  energy  at  frequency  f 

T  wave  period 

Tp  wave  period  corresponding  to  highest  peak  of  the  spectrum 

(Tp  =  f^1) 

e  spectral  width  parameter 

n  sea-surface  elevation 

e  phase 
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ENERGY  SPECTRA  IN  SHALLOW  U.S.  COASTAL  WATERS 

by 

Eduard  F.  Thompson 
I.  INTRODUCTION 

Wind-generated  gravity  waves,  with  periods  ranging  from  1  to  25  sec¬ 
onds,  are  a  primary  concern  in  most  coastal  engineering  work.  Wind  waves 
relate  directly  to  such  processes  as  longshore  sand  movement,  onshore- 
offshore  sand  movement,  and  forces  exerted  on  coastal  structures.  Hence, 
wind  waves  must  be  considered  in  virtually  all  coastal  engineering  en¬ 
deavors,  especially  in  design  and  maintenance  of  harbors  and  exposed 
navigation  channels  and  in  shore  protection  efforts.  Edmisten  (1978) 
presents  an  extensive  list  of  major  applications  for  coastal  wind  wave 
data. 

Coastal  engineers  and  researchers  have  wrestled  for  decades  with  the 
difficult  problem  of  how  to  reduce  the  complexity  of  ocean  wind  waves  to 
simple,  useful  approximations.  The  most  common  approach  has  been  to  rep¬ 
resent  the  dominant  waves  in  a  particular  sea  state  with  three  simple 
parameters:  significant  wave  height,  significant  wave  period,  and  wave 
direction.  The  latter  parameter  is  normally  unavailable,  especially  with 
gage  data;  however,  it  can  be  estimated  from  wind  and  weather  records  for 
days  of  special  interest. 

Practical  techniques  are  available  for  estimating  significant  height 
and  period  from  a  wave  record.  All  of  the  accepted  techniques  produce 
reasonable  estimates  for  significant  height.  However,  estimates  for 
significant  period  are  generally  more  variable  and  their  significance 
is  often  questionable.  One  of  the  best  approaches  to  estimating  sig¬ 
nificant  or  dominant  period  is  to  compute  the  distribution  of  wave  energy 
as  a  function  of  frequency  (frequency  is  the  reciprocal  of  period)  and 
use  the  wave  period  corresponding  to  the  frequency  of  greatest  energy 
density.  This  period  has  a  clear  dynamic  significance.  The  distribution 
of  wave  energy  as  a  function  of  frequency  is  commonly  referred  to  as  the 
wave  energy  epeotrum. 

Although  significant  height,  period,  and  direction  provide  a  useful 
approximation  to  many  sea  states,  these  simple  parameters  omit  details 
of  the  sea  state  which  can  be  important  in  some  applications.  For  ex¬ 
ample,  the  simultaneous  occurrence  of  several  prominent  wave  trains  with 
different  wave  periods  and  directions  is  an  important  consideration  in 
relating  coastal  wave  conditions  to  longshore  sediment  movement.  Design 
of  ships  and  flexible  coastal  structures,  such  as  floating  breakwaters 
and  offshore  platforms,  requires  special  consideration  of  how  much  wave 
energy  can  be  expected  at  frequencies  near  the  resonant  frequencies  of 
the  structure.  Similarly,  effective  harbor  design  requires  that  resonant 
oscillations  in  the  harbor  be  minimized.  Since  the  effectiveness  of 
floating  breakwaters  in  attenuating  wave  energy  depends  on  the  wave  fre¬ 
quency,  floating  breakwater  designs  should  be  formulated  and  evaluated 
with  a  knowledge  of  how  wave  energy  is  distributed  with  frequency. 


The  preceding  examples  illustrate  several  applications  in  which  more 
detailed  specification  of  sea  state  is  needed.  Much  of  this  additional 
needed  detail  is  provided  by  the  energy  spectrum.  Evolving  technology 
and  design  techniques  in  most  areas  of  coastal  engineering  are  generating 
a  similar  need  for  engineers  to  acquire  a  more  advanced  understanding  of 
sea  state  and  its  climate. 

CERC  (and  its  predecessor  the  Beach  Erosion  Board)  has  operated  wave 
gages  at  numerous  U.S.  sites  during  the  last  30  years  and  has  acquired 
a  volume  of  shallow-water  coastal  wave  records  matched  by  few  other  places 
in  the  world.  Most  of  the  gages  were  mounted  on  coastal  piers  in  water 
depths  of  4  to  6  meters.  Significant  wave  heights  and  periods  from  many 
of  the  gages  have  been  summarized  by  Thompson  (1977).  Although  none  of 
the  gages  was  capable  of  sensing  wave  direction,  an  array  of  gages  was 
used  at  one  site  to  obtain  direction  estimates  (Esteva,  1977). 

This  report  presents  spectra  and  spectral  summaries  for  the  selected 
CERC  gage  locations  shown  in  Figure  1;  it  contains  264  high-energy  spec¬ 
tra  which  mostly  represent  shallow-water  wave  conditions.  These  spec¬ 
tra  are  a  significant  addition  to  the  existing  published  literature  on 
shallow-water  spectra.  Field  measurements  continually  show  that  indi¬ 
vidual  spectra  are  quite  erratic  and  variable.  To  identify  general 
spectral  characteristics,  this  report  contains  summaries  of  approxi¬ 
mately  7,000  spectra  and  selected  spectral  parameters;  summaries  of 
parameters  of  the  sea-surface  elevation  distribution  function  are  also 
included. 

The  scope  of  this  report  is  limited.  Fewer  than  one-half  of  the 
CERC  gages  for  which  spectra  are  available  are  included.  Gages  were 
selected  for  optimum  coverage  of  U.S.  coasts  and  optimum  accessibility 
of  archived  spectra.  Only  a  fraction  of  the  available  spectra  from 
each  gage  was  considered.  For  most  gages,  the  spectra  represent  one 
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reasonably  complete  year  of  data.  It  seems  doubtful  that  more  extensive 
summaries  would  be  worth  the  requisite  cost  and  effort  until  techniques 
for  interpreting  and  applying  spectra  are  better  established. 

This  report  is  intended  to  provide  coastal  engineers  a  better  quali¬ 
tative  understanding  of  shallow-water  waves  and  spectra  as  well  as  pro¬ 
viding  quantitative  information  for  engineering  use.  Qualitatively,  it 
provides: 

(a)  Perspective  on  how  to  interpret  shallow-water  wave 
energy  spectra; 

(b)  numerous  examples  of  high-energy,  shallow-water  field 
wave  spectra  from  the  U. S.  Atlantic,  Pacific,  gulf,  and  Great 
Lakes  coasts; 

(c)  evidence  that  the  ratio  of  significant  wave  height  to 
water  depth  never  approaches  the  commonly  used  limit  of  0.78 
in  the  water  depths  considered;  and 

(d)  perspective  on  the  likelihood  of  relatively  high 
individual  waves  occurring  in  succession  in  a  sea  state. 

Quantitatively,  the  report  provides  the  following  information  helpful  to 
coastal  engineers: 

(a)  Average  spectra  corresponding  to  the  more  common  sig¬ 
nificant  height  and  period  combinations  at  each  site; 

(b)  frequency  of  occurrence  of  sea  states  with  major  sec¬ 
ondary  wave  trains; 

(cl  probability  statistics  for  moments  of  the  sea-surface 
elevation  distribution  function,  which  can  be  related  to  prob¬ 
abilities  associated  with  instantaneous  surface  elevations 
above  and  below  the  local  mean;  and 

(d)  a  computer  program  for  identifying  major  peaks  and 
valleys  in  an  irregular  signal;  it  has  many  potential  appli¬ 
cations  including  identification  of  major  peaks  in  a  spectrum 
and  identification  of  meaningful  crests  and  troughs  in  a  time 
series  of  sea-surface  elevation. 

A  brief  discussion  of  the  physical  characteristics  of  shallow-water 
ocean  waves  is  presented  in  Section  II.  Wave  gage  data  collection, 
analysis,  and  summarization  techniques  are  described  in  Section  III 
which  includes  a  discussion  of  the  motivation  for  grouping  spectra  by 
both  significant  height  and  dominant  period  and  computing  average  spec¬ 
tra.  Spectra  computed  for  several  theoretical  cases  of  cnoidal  wave 
profiles  are  presented  for  perspective  in  interpreting  field  spectra. 

The  choice,  computation,  and  interpretation  of  spectral  parameters  and 
sea-surface  elevation  distribution  function  parameters  are  also  dis¬ 
cussed  in  Section  III. 
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The  analyzed  and  summarized  field  wave  spectra  are  discussed  in 
Section  IV  which  includes  samples  of  spectra  representing  the  highest 
measured  wave  conditions  at  each  site,  average  spectra  and  summaries 
of  parameters  derived  from  the  spectrum  and  the  sea-surface  elevation 
distribution  function.  Conclusions  from  the  results  in  Section  IV  are 
provided  in  Section  V.  Major  findings  are  summarized  in  Section  VI. 
Numerous  individual  high-energy  spectral  plots  are  contained  in  Appen¬ 
dix  A.  Average  spectra  are  included  in  Appendixes  B  and  C. 

II.  PHYSICAL  CHARACTERISTICS  OF  SHALLOW- WATER  OCEAN  WAVES 

Wind  waves  change  considerably  when  propagating  from  deep  to  shallow 
water.  Consider  a  train  of  uniform  waves  moving  toward  shore  with  crests 
oriented  parallel  to  the  bottom  contours.  The  wave  height  first  decreases 
slightly  and  then  increases  continually  as  the  waves  move  into  progres¬ 
sively  shallower  water.  The  phase  velocity,  group  velocity,  and  wave¬ 
length  all  decrease  with  decreasing  water  depth.  Wave  energy  is  generally 
reduced  as  waves  propagate  shoreward. 

It  appears  contradictory  that  the  wave  height  in  a  train  of  uniform 
waves  moving  shoreward  increases  while  the  energy  decreases.  Yet  the 
two  effects  occur  together  mainly  because  the  wave  profile  changes.  The 
profile  of  a  wave  in  deep  water  can  usually  be  considered  sinusoidal  un¬ 
less  the  waves  are  exceptionally  high  in  relation  to  their  wavelength. 
However,  the  profile  of  a  steep  wave  in  very  shallow  water  is  decidedly 
nonsinusoidal  (Fig.  2).  The  troughs  are  broad  and  flat,  and  the  crests 
are  narrow  and  high.  Deviation  of  the  wave  profile  from  sinusoidal  be¬ 
comes  important  when  the  relative  water  depth  (ratio  of  water  depth  to 
wavelength)  is  less  than  0.04  for  wind  waves  with  very  low  steepness 
(ratio  of  wave  height  to  wavelength) .  Nonsinusoidal  profiles  also  occur 
for  relative  depths  greater  than  0.04  when  the  wave  steepness  is  greater 
than  about  0.0015;  i.e.,  for  ratios  of  wave  height  to  the  product  of 
gravitational  acceleration  times  wave  period  squared  greater  than  0.00006 
(U.S.  Army,  Corps  of  Engineers,  Coastal  Engineering  Research  Center,  1977, 
p.  2-35).  For  waves  with  very  high  steepness  (greater  than  0.0063  or 
height  divided  by  gT2  greater  than  0.001),  nonsinusoidal  profiles  occur 
in  all  water  depths. 

Since  phase  velocity  of  an  ocean  wave  decreases  with  decreasing  water 
depth,  a  long  wave  crest  approaching  shore  at  an  angle  to  the  bottom  con¬ 
tours  has  a  nonuniforra  phase  velocity  along  the  crest.  The  velocity  dif¬ 
ferences  result  in  an  apparent  bending,  or  refracting,  of  the  crest  which 
decreases  the  angle  between  the  crest  and  the  bottom  contours.  Refraction 
can  affect  a  wave  when  the  relative  water  depth  is  less  than  one-half. 
Relative  depths  for  the  gages  in  this  report  range  typically  from  0.04  to 
1.  Hence,  refraction  has  potentially  affected  many  of  the  wave  measure¬ 
ments  considered  in  this  study.  Refraction  can  act  to  increase  or  de¬ 
crease  wave  heights  depending  on  wave  period  and  direction  and  bottom 
contour  orientation.  Since  most  of  the  gages  are  sufficiently  nearshore 
that  large  waves  can  move  bottom  sediment  seaward  of  the  gage,  the  bottom 
contour  orientation  affecting  refraction  can  change  with  time. 
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Figure  2,  Wave  profiles  and  energy  spectra  for 
several  cnoidal  wave  cases  (record 
length:  512  seconds;  spectral  band¬ 
width:  0.00977  hertz). 

Another  very  important  phenomenon  is  wave  breaking,  which  for  mono¬ 
chromatic  waves  occurs  in  deep  water  approximately  when  the  wave  steep¬ 
ness  exceeds  0.14  or  in  shallow  water  when  the  ratio  of  wave  height  to 
water  depth  exceeds  0.78.  The  heights  of  individual  nearshore  waves  are 
actually  quite  variable  but  generally  conform  to  a  Rayleigh  distribution 
(Thompson,  1974;  Goda  1974;  U.S.  Array,  Corps  of  Engineers,  Coastal  Engi¬ 
neering  Research  Center,  1977).  Thus,  13.5  percent  of  the  waves  can  be 
expected  to  be  higher  than  the  significant  height  and  break  farther  sea¬ 
ward  than  a  wave  with  the  significant  height.  Breaking  of  a  few  very 
large  waves  does  not  visibly  affect  the  significant  height,  but  when  the 
depth-induced  breaking  starts  to  affect  many  large  waves,  the  significant 
height  decreases.  Quantitative  predictions  of  the  decrease  in  significant 
height  as  waves  begin  to  break  nearshore  are  given  by  Seelig  and  Ahrens 
(in  preparation,  1980)  for  a  variety  of  wave  steepnesses  and  beach  slopes. 
Their  design  curves,  based  on  theory  developed  by  Goda  (1975),  indicate 
that  breaking  can  begin  to  affect  significant  height  when  the  ratio  of 
significant  height  to  water  depth  is  greater  than  about  0.5.  This  gen¬ 
eral  observation  is  consistent  with  field  wave  gage  data  presented  by 
Irie  (1975).  For  a  typical  water  depth  at  the  gage  of  5  meters,  break¬ 
ing  might  be  expected  to  have  affected  all  significant  heights  greater 
than  about  2.5  meters. 
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The  large  effect  wave  breaking  can  have  on  significant  height  meas¬ 
ured  at  most  of  the  CERC  gages  is  shown  in  Figure  3.  The  figure  shows 
a  time  history  of  significant  wave  height  for  three  Atlantic  coast  gages 
during  a  large  winter  storm-  The  standard  Nags  Head,  North  Carolina, 
gage  is  represented  along  with  two  staff  gages  on  the  new  CERC  Field 
Research  Facility  (FRF)  pier  30  miles  north  of  Nags  Head.  The  FRF  gage 
farthest  seaward  shows  peak  significant  heights  about  1  meter  higher 
than  the  FRF  gage  390  meters  farther  shoreward.  During  the  peak  of  the 
storm,  significant  heights  at  the  Nags  Head  gage  drop  well  below  heights 
at  the  seaward  FRF  gage.  The  figure  suggests  that  breaking  has  a  major 
influence  on  significant  height  at  both  the  nearshore  FRF  gage  and  the 
Nags  Head  gage. 
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Figure  3.  Time  history  of  significant  wave  height  for  three  North 
Carolina  gages  during  a  large  winter  storm. 


For  some  of  the  gage  sites  included  in  this  report,  and  perhaps  for 
all  except  the  accelerometer  buoy  sites,  a  sandbar  was  often  seaward  of 
the  gage.  The  bar  is  a  somewhat  transient  feature  which  can  have  a  major 
influence  on  the  waves  by  causing  high  waves  to  break  before  they  arrive 
at  the  gages  and  by  acting  to  reflect  energy  seaward  or  trap  it  between 
the  bar  and  shore.  Evidence  presented  by  Wang  and  Yang  (1976)  indicates 
the  bar  acts  mainly  as  an  energy  dissipator  at  low  tide.  At  high  tide 
for  nonbreaking  waves,  the  bar  tends  to  reflect  wave  energy  seaward 
during  onshore  winds  and  to  trap  wave  energy  between  the  bar  and  beach 
during  offshore  winds. 

Since  ocean  waves  in  a  sea  state  are  never  strictly  uniform  in  either 
height  or  frequency,  a  sea  state  is  conveniently  characterized  by  the 
distribution  of  wave  energy  as  a  function  of  frequency  (the  spectrum). 
During  active  wave  growth,  most  energy  from  the  atmosphere  is  added  at 
frequencies  slightly  higher  than  the  dominant  frequency.  The  energy  is 
then  redistributed  by  nonlinear  transfers  within  the  wave  field  to  fre¬ 
quencies  slightly  lower  than  the  dominant  frequency  and  to  very  high 
frequencies  where  it  is  lost  to  breaking  processes.  These  phenomena 
result  in  a  gradual  shift  of  the  spectral  peak  to  lower  frequencies. 
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For  waves  propagating  outside  the  wave-generation  area  the  nonlinear 
energy  transfer  to  low  and  very  high  frequency  continues  without  the 
acconpanying  replenishment  of  energy  at  midfrequencies.  Thus,  swell 
waves  are  characterized  by  narrow  spectra  with  peaks  at  low  frequency. 

There  is  evidence  that  the  intensity  of  nonlinear  energy  transfer 
increases  rapidly  with  decreasing  water  depth.  The  nonlinear  transfer 
may  also  narrow  the  directional  spread  of  energy  in  shallow-water  waves 
(Herterich  and  Hasselmann,  in  preparation,  1980). 

Aerial  photos  of  nearshore  ocean  waters  indicate  remarkable  organi¬ 
zation  in  the  surface  waves  (Fig.  4) ;  the  photos  also  indicate  that  sea 
states  dominated  by  a  single  visible  wave  train  are  relatively  unusual. 
Photos  typically  show  two  or  three  apparently  independent  wave  trains 
with  different,  well-defined  directions  and  wavelengths.  The  relative 
prominence  of  the  trains  often  varies  with  proximity  to  shore  due  to 
frequency-dependent  shoaling  and  refraction,  and  in  some  cases  due  to 
redistribution  of  energy  with  frequency  by  nonlinear  interactions 
between  trains. 

For  example,  although  a  long  swell  wave  train  often  dominates  the 
breaker  zone  in  photos  taken  along  the  southern  California  coast,  the 
train  is  invisible  two  or  three  wavelengths  seaward  of  the  breaker  zone. 
The  swell  waves  are  preferentially  amplified  nearshore.  McClenan  and 
Harris  (1975)  show  several  examples  of  this  process  occurring. 

The  swell  waves  are  also  refracted  more  completely  than  higher  fre¬ 
quency  waves  so  that  at  breaking  they  are  often  nearly  parallel  to  the 
shore.  The  water  particle  velocity  associated  with  low-frequency  swell 
waves  attenuates  slowly  with  depth.  High-frequency  waves  which  have  not 
been  substantially  refracted  can  travel  into  the  swell-created  breaker 
zone  with  a  large  component  of  longshore  motion;  yet,  their  particle 
velocities  attenuate  rapidly  with  depth.  Since  longshore  sediment  move¬ 
ment  requires  both  a  mechanism  for  initiating  sediment  movement  and  a 
longshore  drift  to  transport  it,  two  concurrent  wave  trains  may  move 
sediment  alongshore  more  effectively  than  either  train  acting  alone. 

The  presence  of  a  current  near  a  gage  site  can  alter  the  measured 
characteristics  of  waves.  The  longshore  current  in  and  near  the  breaker 
zone  affects  some  of  the  gage  measurements  in  this  report.  Tidal  cur¬ 
rents  and  wind- generated  currents  may  also  affect  the  measurements.  The 
Gulf  Stream  may  affect  deepwater  waves  approaching  the  U.S.  Atlantic 
coast  by  steepening  waves  enough  to  cause  breaking  or,  in  special  cases, 
by  refracting  the  waves  so  that  they  never  reach  the  shore  (Kenyon,  1971). 
The  rate  of  wave  growth  in  the  vicinity  of  the  current  may  also  be 
affected. 


III.  WAVE  GAGE  DATA  COLLECTION  AND  ANALYSIS 
1.  Collection. 

The  wave  gage  data  used  in  this  report  were  collected  as  part  of  the 
CERC  field  wave  data  collection  program  from  11  gages  at*  the  9  locations 
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shown  in  Figure  1.  Table  1  gives  some  details  of  the  gage  locations, 
gage  types,  and  water  depth  at  the  gages.  The  gage  designs  and  charac¬ 
teristics  are  described  in  detail  by  Thompson  (1977) ,  and  for  the  Great 
Lakes  accelerometer  buoy  installations  in  Thompson  (1978). 

The  continuous -wire  staff  gage  (manufactured  by  the  Baylor  Company, 
Houston,  Texas)  is  the  more  accurate  sensor  because  the  gage  directly 
senses  surface  waves  and  it  interferes  very  little  with  waves  being 
measured.  The  continuous -wire  gage  and  the  pressure  gage  (fabricated 
at  CERC)  give  comparable  surface  wave  height  estimates  when  the  pressure- 
gage  record  is  compensated  for  the  depth-dependent  attenuation  of  the 
dynamic  component  of  pressure  due  to  surface  waves.  Suitably  compen¬ 
sated  pressure-gage  spectra  are  also  comparable  to  continuous -wire  gage 
spectra  except  at  high  frequencies  (Esteva  and  Harris,  1970). 

The  accelerometer-buoy  gage  or  Waverider  (manufactured  by  Datawell, 
Haarlem,  the  Netherlands)  gives  wave  heights  and  spectra  which  are 
reasonably  comparable  to  those  given  by  the  continuous-wire  gage  for 
frequencies  between  0.065  and  0.5  hertz,  including  virtually  all  impor¬ 
tant  frequencies  in  the  Great  Lakes  (Pitt,  Driver,  and  Ewing,  1978). 

The  step-resistance  staff  gage  (fabricated  at  CERC)  has  been  shown  to 
have  a  consistent  bias  toward  hi$i  wave  heights  (Esteva  and  Harris, 

1970) .  Significant  height  estimates  from  the  step-resistance  gage  may 
be  30  centimeters  too  high  during  low  wave  conditions  and  20  percent  too 
high  during  high  wave  conditions.  Aside  from  the  tendency  for  over¬ 
estimating  wave  height  and  energy,  Esteva  and  Harris  show  evidence  that 
spectra  for  step-resistance  gage  records  are  reasonably  consistent  with 
continuous-wire  gage  spectra. 

Time  series  of  sea-surface  elevation  from  staff  gages,  subsurface 
pressure  from  pressure  gages,  and  double-integrated  sea-surface  vertical 
acceleration  from  buoy  gages  were  transmitted  by  telephone  line  to  the 
CERC  laboratory  and  used  to  generate  pen-and-ink,  strip-chart  records 
and  digital  records  on  computer-compatible,  seven-track  magnetic  tape. 
Digital  records  are  20  minutes  long.  A  digital  record  from  a  single 
station  contains  data  points  at  0.25-second  intervals  giving  a  total  of 
4,800  data  points  per  location  per  20-minute  record.  Further  details  of 
CERC  field  wave  data  collection  procedures  are  given  in  Thompson  (1974, 
1977) . 

2.  Analysis. 

a.  Spectral  Analysis.  The  digital  field  wave  gage  records  were 
analyzed  by  the  CERC  routine  analysis  system  (Thompson,  1977).  One 
record  per  6  hours  is  routinely  analyzed;  however,  one  record  per  2 
hours  was  analyzed  for  some  of  the  Great  Lakes  data.  The  routine  anal¬ 
ysis  times  are  chosen  to  approximately  coincide  with  synoptic  times 
(0100,  0700,  1300,  and  1900  e.s.t.).  The  procedure  includes  editing 
each  20-minute  time  series  to  eliminate  data  points  which  are  obviously 
bad,  computing  and  testing  the  distribution  function  of  the  edited  data 
points,  and  computing  the  variance  or  energy  spectrum  for  the  time  series 
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smoothed  with  a  cosine  bell  data  window.  The  spectrum  is  computed  with 
a  fast  Fourier  transform  (FFT)  algorithm  applied  to  4,096  data  points 
(17-minute  and  4-second  record). 

The  spectral  analysis  procedure  assigns  a  fraction  of  the  total  var¬ 
iance  of  the  record  to  each  of  1,024  frequencies,  or  frequency  lines. 

It  is  difficult  to  deal  with  such  finely  resolved  spectra.  To  gain  sta¬ 
tistical  stability,  groups  of  11  successive  lines  in  each  spectrum  are 
combined  into  bands.  The  resultant  band  width  is  0.01074  hertz  over  the 
full  range  of  frequencies  considered  (0.03  to  1.00  hertz).  The  uniform 
spectral  resolution  over  the  frequency  range  gives  a  nonuniform  resolu¬ 
tion  in  wave  period.  Significant  height  is  estimated  as  four  times  the 
square  root  of  the  total  variance  assigned  to  wave  frequencies  of  inter¬ 
est,  usually  0.03  to  1.00  hertz. 

The  frequency-dependent  responses  of  some  of  the  gage  types  required 
some  special  treatment  of  the  spectra.  Accelerometer-buoy  gage  spectra 
were  terminated  at  the  high-frequency  end  at  0.S  hertz.  A  special  low- 
frequency  bound  of  about  0.065  hertz  on  the  spectrum  would  have  been 
appropriate  but  was  not  used  because  energy  at  frequencies  between  that 
and  the  low-frequency  cutoff  for  other  gage  types  (0.03  hertz)  was  in¬ 
significant  in  the  Great  Lakes  buoy  records. 

Pressure-gage  records  were  compensated  for  hydrodynamic  attenuation 
of  the  pressure  signal  with  depth  by  multiplying  the  energy  in  each 
spectral  band  by  a  frequency-dependent  factor.  The  compensation  factor 
becomes  very  large  at  high  frequencies  where  the  small  pressure  signal 
is  often  obscured  by  noise.  Hence,  it  is  necessary  to  terminate  the 
high-frequency  end  of  pressure-gage  spectra  at  a  frequency  below  the 
standard  1-hertz  cutoff  for  staff  gages.  The  high-frequency  cutoff 
used  for  the  Michigan  City  and  Presque  Isle  pressure  gages  is  0.33 
hertz.  The  cutoff  for  the  Pt.  Mugu  pressure  gage  is  0.31  hertz.  All 
spectral  energy  assigned  to  frequencies  above  the  cutoffs  was  omitted 
from  both  the  spectrum  and  the  significant  height  estimate.  The  neg¬ 
lected  energy  is  unimportant  during  high  wave  conditions,  but  it  can 
increase  significant  height  by  as  much  as  0.5  meter  during  low  recorded 
wave  conditions.  The  latter  effect  is  often  partially  balanced  by  slight 
overcompensation  of  the  high-frequency  energy  retained  in  the  spectrum. 

b.  Spectral  Summarization. 

(1)  Averaging.  Individual  shallow-water  ocean  wave  energy 
spectra  are  quite  irregular  and  can  change  substantially  and  somewhat 
erratically  with  time  (see  Fig.  5).  Because  of  these  problems,  it  is 
difficult  to  identify  typical  spectra  for  each  wave  gage  site.  However, 
it  is  desirable  to  condense  the  masses  of  individual  spectra  into  a  more 
concise  form. 

At  this  time,  no  completely  satisfactory  technique  is  in  use  for 
isolating  and  summarizing  general  characteristics  of  field  spectra.  One 
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Figure  5.  Sequence  of  spectra  computed  at  6-hour  intervals  for 
Nags  Head,  North  Carolina,  beginning  at  0642  e.s.t., 
20  April  1969.  Hs  is  190  centimeters  for  the  top 
spectrum  decreasing  to  97  centimeters  for  the  bottom 
spectrum.  Major  spectral  peaks  are  marked  with 
asterisks. 


common  approach  is  to  compute  spectral  statistics  which  are  often  based 
on  moments  of  the  spectrum,  where  the  nth  spectral  moment,  is 

defined  as 


%  =  /  S(f)  f*df  (1) 

o 

where  f  is  the  frequency  and  S(f)  the  spectral  energy  at  frequency  f 
This  approach  is  used  in  many  cases  because  it  is  the  only  tractable  way 
to  deal  with  large  numbers  of  spectra.  However,  some  of  these  spectral 
statistical  parameters  have  disadvantages  which  will  be  discussed  later. 

A  relatively  new  and  promising  approach  to  spectral  parameterization 
involves  the  use  of  eigenfunctions  (Vincent  and  Resio,  1977). 

Another  approach  to  the  problem  of  how  to  characterize  and  summarize 
ocean  wave  spectra  has  been  used  in  conjunction  with  ship-response  pre¬ 
diction.  Hoffman  (1974)  organized  deepwater  spectra  into  10  groups 
from  a  site  in  the  North  Atlantic  according  to  significant  height.  For 
each  group,  the  mean  and  standard  deviation  of  spectral  energy  in  each 
band  were  computed.  The  average  energy  of  the  one-third  highest  and 
one-third  lowest  values  in  each  band  was  also  computed.  The  importance 
of  accounting  for  variability  of  spectral  shape  when  predicting  ship¬ 
bending  moment  responses  was  demonstrated. 

An  averaging  procedure  for  summarizing  deepwater  spectra  from  a  site 
in  the  North  Atlantic  was  also  used  by  Gospodnetic  and  Miles  (1974). 
Spectra  were  grouped  according  to  both  significant  wave  height  and 
average  wave  period.  The  average  spectra  were  made  dimensionless  to 
facilitate  comparison  of  spectral  shapes.  Although  most  of  the  varia¬ 
tions  were  removed  by  using  the  dimensionless  format,  it  was  concluded 
that  even  the  dimensionless  spectra  vary  systemmatical ly  with  both  height 
and  period. 

In  very  shallow  water,  spectral  characteristics  would  be  expected  to 
have  a  pronounced  systematic  relationship  to  both  significant  height  and 
dominant  spectral  period,  althougn  the  relationship  would  be  somewhat 
obscured  when  major  secondary  wave  trains  exist.  Wave  height  is  impor¬ 
tant  because  high  waves  lead  to  wave  breaking  in  shallow  water.  Thus, 
any  reasonable  attempt  to  compute  average  spectra  at  a  shallow-water 
location  must  include  stratification  of  the  spectra  according  to  wave 
height  or  energy. 

The  extent  of  wave  shoaling  is  specified  by  the  relative  water  depth, 
d/Lp,  which  is  directly  related  to  the  dominant  wave  period.  Waves  in 
very  shallow  water  assume  a  nonsinusoidal  profile  which  gives  rise  to 
spectra  with  peaks  at  multiples  of  the  dominant  frequency.  An  example 
of  the  profile  and  spectrum  for  several  cnoidal  wave  cases  is  shown  in 
Figure  2;  another  example  is  an  aerial  photo  spectrum  pair  in  McClenan 
and  Harris  (1975,  p.  63-64).  Because  of  such  systematic  effects  on  the 
shallow-water  spectrum,  spectra  should  also  be  stratified  by  d/I p  or 
by  a  characteristic  wave  period. 
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Deviations  of  the  wave  profile  from  sinusoidal  occur  even  in  deep 
water  for  very  steep  waves  as  discussed  previously.  Spectra  for  such 
waves  can  also  be  expected  to  show  energy  at  multiples  of  the  dominant 
frequency.  In  such  cases  the  spectral  form  is  related  to  both  wave 
height  and  period. 

Because  of  expected  systematic  influences  of  wave  height  and  period 
on  the  shallow-water  spectrum,  spectra  in  this  study  were  grouped  into 
sets  by  both  significant  height  and  period  corresponding  to  the  highest 
spectral  peak.  Significant  height  and  peak  period  are  convenient  because 
they  are  comparable  to  height  and  period  parameters  used  extensively  in 
past  nonspectral  wave  analyses. 

Despite  the  grouping,  major  differences  between  individual  spectral 
shapes  in  each  set  can  still  occur.  Differences  between  spectra  at  one 
location  for  several  cases  in  which  the  significant  height  and  peak 
period  were  near  the  annual  means  are  shown  in  Figure  6.  The  differences 
in  spectral  shape  between  individual  spectra  in  a  set  were  greater  than 
the  differences  in  the  mean  spectrum  between  sets  in  this  study  and  in 
Hoffman's  (1974)  study.  However,  any  systematic  characteristics  of  each 
set  should  emerge  in  the  mean  spectrum  and  standard  deviation  if  the 
number  of  samples  is  sufficient.  This  does  not  imply  that  the  mean 
spectra  necessarily  represent  characteristics. 


Figure  6.  Several  individual  spectra  in  the  same  Hs  -  T p  group 

(Hg  =  91  to  122  centimeters;  T p  =  8.5  to  9.3  seconds). 
Nags  Head,  North  Carolina. 


Spectra  were  grouped  by  30-centimeter  intervals  of  significant 
height  and  by  variable  intervals  of  peak  period  (Table  2).  Variable 
period  intervals  were  needed  because  energy  is  computed  for  nonuniform 
period  bands  (but  uniform  frequency  bands)  in  the  CERC  spectral  analysis 
program.  Selection  of  longer  period  intervals  was  necessary  so  that 
each  interval  encompasses  only  one  spectral  band.  If  the  intervals  had 
included  more  than  one  spectral  band,  the  averaging  procedure  would  tend 
to  produce  an  unrepresentatively  low  and  broad  peak  in  the  average 
spectrum. 


Table  2.  Peak  period  grouping  intervals. 


Period  interval 

Frequency  bands 

Period  interval 

Frequency  bands 

in  interval 

in  interval 

(s) 

(No.) 

(s) 

(No.) 

1  to  2 

47 

7.8  to  8.5 

1 

2  to  3 

IS 

8.5  to  9.3 

1 

3  to  4 

8 

9.3  to  10.3 

1 

4  to  5 

4 

10.3  to  11.6 

1 

5  to  6 

4 

11.6  to  13.3 

1 

6  to  6*6 

1 

13.3  to  15.5 

1 

6.6  to  7.2 

1 

15.5  to  18.6 

1 

7.2  to  7.8 

1 

18.6  to  23.3 

1 

Period  intervals  of  1  second  were  used  for  periods  shorter  than  6 
seconds  (Table  2).  Individual  spectra  with  relatively  short  peak  pe¬ 
riods  are  generally  more  poorly  focused  in  frequency  than  spectra  with 
long  peak  periods.  Thus,  period  intervals  encompassing  more  than  one 
spectral  band  could  be  used  for  short  periods  without  creating  much 
artificial  broadening  of  the  peak  in  the  average  spectrum.  This  treat¬ 
ment  is  further  defensible  on  the  grounds  that  (a)  short -period  waves 
are  in  relatively  deeper  water  than  long-period  waves  at  the  measure¬ 
ment  site  and  hence  exhibit  less  depth- induced  modification,  and  (b) 
spectra  with  peak  periods  shorter  than  6  seconds  are  generally  not  im¬ 
portant  for  coastal  engineering  design  applications  at  ocean  sites. 

After  the  spectra  were  grouped  into  height-period  intervals,  average 
spectra  were  confuted.  Intervals  were  selected  for  averaging  to  include 
all  intervals  containing  a  reasonable  number  of  cases  and  all  high  wave 
intervals  regardless  of  the  number  of  cases.  The  minimum  number  of  cases 
considered  for  low  and  moderate  waves  was  10  at  most  locations.  The  wave 
height  above  which  all  intervals  were  considered  varies  between  locations. 
The  average  spectrum  for  a  height-period  interval  was  obtained  by  aver¬ 
aging  energy  values  in  each  band.  The  standard  deviation  of  spectral 
energy  values  about  the  mean  for  each  band  was  also  computed.  Averages 
and  standard  deviations  of  spectra  for  the  locations  listed  in  Table  1 
are  discussed  in  the  next  section. 

It  is  important  to  note  that  the  spectrum  for  any  individual  record 
never  duplicates  the  mean  spectrum.  Any  application  which  is  sensitive 
to  the  precise  form  of  the  spectrum  should  make  use  of  individual  spectra 
rather  than  the  mean  spectrum.  One  reasonable  approach  might  be  to  select 
a  small  random  sample  of  spectra  from  each  group,  apply  each  spectrum  to 
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the  problem  at  hand,  and  combine  the  individual  results  to  estimate 
response  probabilities,  as  suggested  by  Hoffman  (1974,  1975). 

(2)  Spectral  Parameters. 

(a)  Number  of  Peaks.  It  is  relatively  common  in  the  ocean 
for  two  or  more  independent  wave  trains  with  different  frequency  and 
direction  to  occur  simultaneously  (e.g.,  Harris,  1972;  McClenan  and 
Harris,  1975;  Ochi  and  Hubble,  1976;  Thompson,  1977).  The  marine  sur¬ 
face  observation  reporting  code  even  includes  provisions  for  reporting 
major  secondary  wave  trains  observed  visually.  Certainly,  the  concept 
of  a  simple,  single-peaked  spectrum  is  misleading  in  many  ocean  wave 
records.  Yet,  there  is  a  lack  of  reliable  quantitative  information  as 
to  how  often  major  secondary  trains  occur  and  how  much  energy  they 
contain. 

To  investigate  the  occurrence  of  multiple  wave  trains,  a  computer 
routine  was  adapted  for  use  in  identifying  major  spectral  peaks.  The 
routine  (described  in  App.  D)  smoothes  over  spectral  peaks  and  valleys 
which  differ  in  energy  density  by  less  than  3  percent  of  the  total 
energy  in  the  spectrum.  Examples  of  individual  spectra  were  shown  in 
Figure  5.  The  major  peaks  identified  in  each  case  by  the  computer  rou¬ 
tine  are  marked  with  asterisks. 

A  fraction  of  the  spectral  energy  is  assigned  to  each  major  spectral 
peak  by  a  simple  method.  The  spectrum  is  partitioned  at  the  lowest  point 
between  successive  major  peaks  (Fig.  7).  All  spectral  energy  within  the 
partition  for  a  peak  is  assigned  to  that  peak.  The  energy  is  expressed 
as  a  dimensional  quantity  and  as  a  fraction  of  the  total  spectral  energy. 


Figure  7.  Technique  for  partitioning  spectral  energy 
and  assigning  the  energy  to  major  peaks. 
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Statistics  on  the  occurrence  of  multiple  spectral  peaks  identified 
by  the  above  techniques  (presented  in  the  next  section)  are  based  only 
on  the  part  of  the  spectrum  between  0.03  and  0.5  hertz  for  the  staff  and 
accelerometer-buoy  gages.  The  high-frequency  cutoff  for  the  pressure 
gages  was  0.33  hertz  in  the  Great  Lakes  and  0.31  hertz  at  Pt.  Mugu  as 
discussed  earlier. 

(b)  Spectral  Peakedness.  It  is  highly  desirable  to  param¬ 
eterize  a  spectrum.  Various  parameters  have  been  considered  in  the  lit¬ 
erature,  but  most  are  derived  from  moments  of  the  spectrum  as  defined  in 
equation  (1). 

In  practical  application  equation  (1)  becomes 


rHF 

m n  =  J  S(f)  f"df  (2) 

fLF 

where  f pp  is  the  low-frequency  spectral  cutoff,  and  fyp  the  high- 
frequency  spectral  cutoff  used  for  computation.  When  applied  to  field 
records  of  finite  length,  the  integral  in  equation  (2)  becomes  a 
summation. 

Equation  (2)  shows  that  the  higher  the  value  of  fpp,  the  greater 
the  magnitude  of  provided  that  S(f)  is  greater  than  zero  at  high 

frequencies.  In  practice,  the  choice  of  fpp  has  a  significant  influ¬ 
ence  on  the  values  computed  for  m n,  especially  for  the  higher  moments. 
The  choice  of  fpp  is  less  ambiguous  in  practice  because  contributions 
to  ra^  from  the  low-energy  and  low-frequency  part  of  the  spectrum  are 
generally  negligible. 

The  sensitivity  of  m„  and  parameters  based  on  to  the  choice 

of  flip  was  investigated  by  Rye  (1977).  He  considered  two  theoretical 
spectral  shapes:  a  sharply  peaked  JONSWAP  spectrum  and  a  flat  JONSWAP 
spectrum  equivalent  to  a  Pierson-Moskowitz  spectrum.  Rye  recommended 
three  spectral  parameters  for  general  use: 

(a)  Significant  wave  height,  Hg  where 

Hs  =  (3) 

(b)  period  of  the  spectral  peak,  T 'p 

(c)  spectral -peakedness  parameter,  Qp  where 

00 

Qp  =  ^  /f[s(f)]2df  (4) 

The  spectral-peakedness  parameter  was  originally  proposed  and  shown  to 
be  directly  related  to  the  average  number  of  high  wavesj  oocurring  in 


27 


succession  in  simulated  wave  records  by  Goda  (1970).  However,  Goda  (1976) 
could  not  find  a  clear  relationship  between  Qp  and  the  extent  of  group¬ 
ing  of  high  waves  in  field  records. 

The  value  of  Qp  as  a  spectral  parameter  is  shown  in  Figure  8.  The 
two  spectra  shown  have  nearly  the  same  significant  height  and  peak  period. 
The  spectra  also  have  comparable  values  of  the  frequently  used  spectral- 
width  parameter,  e,  defined  as 


=  1  - 


However,  Qp  for  one  spectrum  is  66  percent  greater  than  for  the  other, 
indicating  substantial  differences  in  spectral  shape. 


On  =3.1 


Gn=f.9 


Figure  8. 


Frequency 

Comparison  of  two  measured  spectra 
from  the  North  Atlantic  Ocean  with 
He  =3.3  meters,  Tp  =  10.5  seconds, 
and  e  =  0.61  (after  Hoffman,  1974). 


Both  Qp  and  e  are  designed  to  fit  the  concept  of  a  single-peaked 
spectrum,  but  multipeaked  ocean  wave  spectra  are  commonly  observed.  The 
suitability  of  spectral  parameters  for  representing  multipeaked  field 
wave  spectra  has  received  little  attention  in  the  literature,  although 
Ochi  and  Hubble  (1976)  have  developed  a  computationally  complicated  tech¬ 
nique  for  parameterizing  field  spectra  with  two  major  peaks. 

Since  Qp  seems  to  be  a  potentially  important  and  useful  parameter, 
it  was  computed  for  each  of  the  spectra  considered  in  this  report.  Qp 
was  based  on  the  part  of  the  spectrum  between  0.03  and  0.2  hertz  for 
staff  and  accelerometer-buoy  gages,  between  0.03  and  0.33  hertz  for  the 
Great  Lakes  pressure  gages,  and  between  0.03  and  0.31  hertz  for  the  Pt. 
Mugu  pressure  gage. 

c.  Parameters  of  Distribution  Function  of  Sea-Surface  Elevation. 

The  distribution  function  for  instantaneous  sea-surface  elevations  pro¬ 
vides  useful  insight  on  shallow-water  spectra  and  spectral  parameters. 

It  also  provides  probabilities  associated  with  instantaneous  surface 
elevations  above  the  mean. 

The  distribution  function  of  sea-surface  elevations  can  conveniently 
be  parameterized  by  its  moments  defined  as: 

%  '  jj "?  <« 

where 

<\n  =  nth  moment  of  the  distribution  function  of  sea-surface 

elevations 

N  =  number  of  intervals  in  the  distribution  function 

r\^  =  sea-surface  elevation  associated  with  the  ith  interval 

in  the  distribution  function 

p(n  )  =  probability  associated  with  n,* 

i  ^ 

The  zeroth  and  first  moments,  qQ  and  qj ,  are  equivalent  to  the  mean 
and  variance  of  the  distribution  function.  q3  and  q4  are  often  re¬ 
ferred  to  as  the  skewness  and  kurtosis  of  the  distribution  function. 

The  distribution  function  of  sea-surface  elevations  is  often  assumed 
to  be  Gaussian.  When  normalized,  the  Gaussian  distribution  function  has 
a  mean  of  zero,  a  variance  of  one,  a  skewness  of  zero,  and  a  kurtosis  of 
three . 

Steep  waves  in  shallow  water  assume  a  decidedly  nonsinusoidal  pro¬ 
file  with  broad,  flat  troughs  and  narrow,  high  crests.  The  distribution 
of  sea-surface  elevations  measured  at  a  point  is  obviously  non-Gaussian. 
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The  skewness  and  perhaps  the  kurtosis  of  surface  elevations  in  very 
shallow  water  can  be  expected  to  differ  from  the  values  for  a  normal¬ 
ized  Gaussian  distribution.  For  the  cnoidal  wave  profile  in  Figure 
9,  the  skewness  and  kurtosis  are  equal  to  2.0  and  5.8,  respectively. 

The  figure  also  shows  the  distribution  function,  skewness,  and  kurtosis 
for  several  idealized  wave  profiles. 


Wove  Profile 


Distribution  Function  of 
Seo-Surfoce  Elevotions 


Figure  9. 


Distribution  function  of  sea- 
surface  elevations  for  several 
idealized  wave  profiles. 


Both  skewness  and  kurtosis  for  the  normalized  distribution  function 
of  sea-surface  elevations  were  computed  for  each  record  and  are  summa¬ 
rized  in  the  following  section.  For  the  accelerometer-buoy  data,  the 
distribution  functions  represent  doubly  integrated  vertical-surface 
accelerations  rather  than  directly  measured  sea-surface  elevations. 

For  the  pressure-gage  data,  the  distribution  functions  represent  sub¬ 
surface  pressure  fluctuations. 


IV.  SPECTRA,  SPECTRAL  SUMMARIES,  AND  SEA-SURFACE 
ELEVATION  DISTRIBUTION  FUNCTION  PARAMETERS 

1.  High-Energy  Wave  Spectra. 

Since  high-energy  spectra  are  of  greatest  concern  in  most  engineering 
applications,  plots  of  individual  spectra  containing  the  most  energy  for 
each  station  are  provided  in  Appendix  A.  Spectra  are  grouped  by  gage 
location  in  the  order  listed  in  Table  1.  The  appendix  includes  24  spec¬ 
tra  for  each  location,  arranged  and  numbered  in  descending  order  of  sig¬ 
nificant  height. 

Each  plot  in  Appendix  A  shows  the  significant  wave  height,  period 
corresponding  to  the  highest  spectral  peak,  spectral-peakedness  parameter 
((^>),  and  number  of  major  spectral  peaks  identified  by  the  computer  rou¬ 
tine  SMOOTH  (App.  D).  Spectral  plots  for  days  associated  with  tropical 
storms  and  hurricanes  are  labeled  accordingly.  The  energy  density  scale 
can  vary  between  plots. 

Pressure-gage  spectra  in  Appendix  A  are  cut  off  at  high  frequencies 
above  0.33  hertz  for  the  Great  Lakes  gages  and  0.31  hertz  for  the  Pt. 

Mugu  gage.  Several  of  the  Pt.  Mugu  high-energy  spectra  were  obviously 
overcompensated  for  attenuation  of  pressure  with  depth  so  that  the  high- 
frequency  energy  was  much  greater  than  the  low-frequency  energy.  These 
spectra  were  omitted  from  Appendix  A.  This  problem  can  arise  because 
the  high-frequency  cutoff  is  constant  while  the  compensation  increases 
as  the  depth  of  water  over  the  gage  increases  with  the  tide.  Thus,  the 
high-frequency  end  of  the  spectrum  is  sometimes  overcompensated  when 
water  levels  are  high.  This  problem  did  not  arise  with  pressure-gage 
spectra  in  the  Great  Lakes  (where  there  is  no  tide) . 

Energy  computed  for  wave  periods  longer  than  30  seconds  (frequencies 
less  than  0.043  hertz)  does  not  directly  represent  wind- generated  waves. 

In  some  cases  it  is  questionable  whether  such  very  long-period  energy 
has  any  physical  meaning.  Therefore,  it  was  omitted  from  spectral  plots 
for  the  Great  Lakes  and  Pacific  coast  locations  and  should  be  ignored  in 
plots  for  the  Atlantic  and  gulf  coast  locations. 

Spectral  plots  for  each  location  in  Appendix  A  are  followed  by  a  table 
which  contains  auxiliary  data  for  each  plot,  including  sequential  number, 
date,  time,  significant  wave  height,  peak  wave  period,  water  depth,  rela¬ 
tive  water  depth  (where  Lp  is  the  shallow-water  wavelength  for  a  wave 
with  frequency  equal  to  the  frequency  of  maximum  spectral  energy  density) , 
ratio  of  significant  wave  height  to  water  depth,  wave  steepness,  ratio  of 
depth  to  gT$  (which  is  related  to  wavelength),  spectral-peakedness  pa¬ 
rameter,  ana  skewness  and  kurtosis  of  the  distribution  function  of  sea- 
surface  elevations.  Wavelength  used  in  the  tables  was  computed  in  each 
case  from  period  corresponding  to  the  highest ' spectral  peak  and  the  equa¬ 
tions  of  linear  wave  theory.  The  largest  ratio  of  significant  wave  height 
to  depth  is  0.54.  None  of  the  significant  height-to-depth  ratios  approach 
the  commonly  used  limit  of  0.78. 
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Multiple  major  spectral  peaks  are  evident  in  many  of  the  spectra  in 
Appendix  A.  However,  some  of  the  multipeaked  spectra  and  many  of  the 
single-peaked  spectra  appear  to  approximate  the  form  of  the  JONSWAP  spec¬ 
trum  for  waves  undergoing  active  generation  (Hasselmann,  et  al.,  1973). 
Only  seven  of  the  spectra  selected  for  inclusion  in  Appendix  A  were  col¬ 
lected  during  hurricanes.  This  small  sample  does  not  illustrate  any 
clear  differences  between  spectra  generated  by  hurricanes  and  spectra 
generated  by  extratropical  storms. 

2.  Mean  Spectra. 

a.  Grouped  by  Height  and  Period.  Spectra  from  each  gage  site  were 
grouped  into  sets  by  significant  height  and  peak  spectral  period.  Inter¬ 
vals  of  30.5  centimeters  were  used  for  significant  height;  intervals  of 
about  1  second  were  used  for  peak  period  (Table  2). 

Mean  spectra  were  plotted  for  all  height-period  intervals  containing 
more  than  one  case  with  heights  greater  than  a  cutoff  height,  which 
varied  with  location.  In  addition,  mean  spectral  plots  were  generated 
for  height-period  intervals  with  heights  below  the  cutoff  if  there  were 
a  reasonably  large  number  of  cases  in  the  interval.  This  approach  pro¬ 
vided  mean  spectral  plots  for  the  most  statistically  meaningful  low  wave 
cases  and  for  all  high  wave  cases.  Energy  computed  for  periods  longer 
than  30  seconds  was  omitted  from  plots  for  Great  Lakes  and  Pacific  coast 
locations  as  discussed  earlier. 

The  mean  spectral  plots  grouped  by  location  are  given  in  Appendix  B. 
The  approximate  relative  water  depth  and  ratio  of  significant  wave  height 
to  water  depth  for  each  location  are  listed  in  Table  B-l;  the  dates 
covered  are  in  Table  1. 

The  range  of  significant  wave  height  and  peak  spectral  period  repre¬ 
sented  and  the  number  of  cases  are  listed  on  each  plot.  Each  plot  repre¬ 
senting  more  than  one  case  contains  several  dashlines  in  addition  to  the 
solid  line  showing  the  mean  spectrum.  The  dashline  nearest  the  mean  and 
above  it  represents  the  mean  plus  one  standard  deviation.  Similarly, 
the  dashline  nearest  the  mean  and  below  it  represents  the  mean  minus  one 
standard  deviation.  When  the  mean  minus  one  standard  deviation  is  nega¬ 
tive,  it  is  plotted  as  zero.  For  plots  representing  more  than  two  cases, 
another  dashline  above  the  mean  shows  the  highest  single-energy  density 
value  occurring  in  each  band.  This  line  forms  an  envelope  inside  which 
all  of  the  individual  spectra  would  fall.  The  dashlines  should  not  be 
treated  as  a  spectrum,  and  are  included  only  as  indicators  of  the  observed 
variability  of  energy  density  in  each  spectral  band.  Some  of  the  plots 
show  a  relatively  large  concentration  of  energy  at  the  lowest  frequency 
(0.005  hertz).  This  concentration  is  at  a  frequency  which  is  irrelevant 
to  wind- generated  gravity  waves  and  should  be  ignored. 

The  mean  spectral  plots  for  each  location  in  Appendix  B  are  arranged 
in  order  of  increasing  significant  height  and  peak  period  in  most  cases; 
i.e.,  all  mean  spectra  in  a  column  represent  the  same  peak-period  interval 
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and  all  mean  spectra  in  a  row  represent  the  same  significant  height 
interval.  Thus,  the  effect  of  increasing  height  on  spectra  with  similar 
peak  periods  and  the  effect  of  increasing  peak  period  on  spectra  with 
similar  heights  can  be  easily  followed.  A  small  number  of  plotted  mean 
spectra  are  not  included  in  Appendix  B  because  the  spectra  (a)  usually 
represented  low  or  moderate  wave  heights  with  relatively  few  cases  or 
high  wave  heights  with  only  one  case  for  which  the  spectrum  is  shown  in 
Appendix  A,  or  (b)  were  high-energy  cases  for  Pt.  Mugu  which  had  obviously 
been  overcoaipensated  for  attenuation  of  pressure  with  depth. 

Mean  spectra  for  low  wave  heights  and  peak  periods  less  than  5  seconds 
along  the  Atlantic  coast  show  evidence  of  secondary  peaks  at  low  frequency 
in  many  cases.  The  low-frequency  peaks  are  broad  and  flat,  are  generally 
located  near  0.1  hertz,  and  are  likely  to  represent  low-energy  swell  in 
the  Atlantic  Ocean.  Corresponding  broad,  low-frequency  secondary  peaks 
are  also  evident  in  mean  spectra  for  Pt.  Mugu  on  the  Pacific  coast.  These 
peaks  are  not  evident  in  mean  spectra  for  the  gulf  and  Great  Lakes  gages, 
although  the  Presque  Isle  pressure-gage  spectra  show  a  sharp  secondary 
peak  which  is  consistently  centered  on  0.038  hertz.  This  low-frequency 
peak  at  Presque  Isle  may  not  be  physically  meaningful.  For  all  locations, 
the  differences  between  mean  spectra  are  smaller  than  the  differences 
between  the  individual  spectra  used  to  compute  the  means. 

Mean  spectra  from  the  ocean  gages  for  peak  periods  longer  than  about 
10  seconds  show  evidence  of  secondary  peaks  at  high  frequency,  especially 
for  high  wave  heights.  For  example,  several  selected  mean  spectra  for 
high  wave,  long-period  groups  are  plotted  in  Figure  10  in  dimensionless 
form.  Energy  density  was  made  dimensionless  for  each  mean  spectrum  by 
dividing  by  the  total  energy  in  the  spectrum  and  converting  to  percent; 
frequency  was  made  dimensionless  by  dividing  by  the  peak  frequency. 

The  figure  clearly  shows  that  the  main  secondary  peak  in  each  mean 
spectrum  is  at  twice  the  peak  frequency.  The  spectrum  for  one  cnoidal 
wave  (see  Fig.  2)  is  also  shown  to  clearly  resemble  the  field  spectra. 

A  section  of  the  I.ake  Worth  pen-and-ink  strip-chart  record  taken  at 
nearly  the  same  time  as  the  Lake  Worth  spectrum  in  Figure  10  is  shown 
in  Figure  11.  The  profiles  of  the  larger  waves  closely  resemble  the 
cnoidal  wave  profile.  Many  of  the  high-frequency  secondary  peaks  in 
mean  spectra  for  high,  long-period  waves  are  probably  indicative  of 
cnoidal-type  wave  profiles.  Spectra  for  high,  long-period  waves  in 
which  secondary  peaks  at  about  two  or  three  times  the  peak  frequency 
represent  independent  wave  trains  are  expected  to  be  rare. 

The  plots  in  Appendix  B  also  show  interesting  variations  with  wave 
height  in  mean  spectra  with  short  peak  periods.  In  some  cases,  the 
mean  spectra  for  high  waves  exhibit  secondary  peaks  at  low  frequency. 

The  emergence  of  the  low-frequency  secondary  peaks  for  Nags  Head  is 
shown  by  dimensionless  spectra  plotted  in  Figure  12.  The  figure  also 
indicates  that  the  largest  secondary  peak  is  at  one-half  the  frequency 
of  the  dominant  frequency.  A  section  of  pen-and-ink  strip-chart  record 
corresponding  to  two  high  wave,  short-period  cases  is  shown  in  Figure  13. 
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Figure  12.  Dimensionless  near  spectra  for  several  high  wave, 
short -period  grojps  at  Nags  Head,  North  Carolina. 
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Figure  13.  Two  pen-and-ink  strip-chart  records  and  corresponding  spectra  frosi 
Nags  Head,  North  Carolina,  (a)  Record  at  0800  e.s.t.,  20  April  1969 
analysis  at  0642  e.s.t.,  20  April  1969;  (b)  record  at  2000  e.s.t., 

1  March  1969;  analysis  at  2040  e.s.t.,  1  March  1969. 


In  one  case  the  spectrum  has  a  single  dominant  peak;  in  the  other  case 
the  spectrum  has  a  large  peak  at  one-half  the  dominant  frequency  and 
another  major  peak  at  about  1.5  times  the  dominant  frequency. 

A  laboratory  parallel  of  this  behavior  has  been  observed  by  Sawaragi 
and  Iwata  (1976).  A  single-peaked  spectrum  was  generated  in  shallow 
water  of  uniform  depth.  As  the  waves  propagated  and  broke,  an  increas¬ 
ing  fraction  of  the  energy  shifted  to  frequencies  which  are  multiples  of 
the  dominant  frequency,  fp.  Eventually,  the  spectral  peak  at  2fp  was 
higher  than  the  peak  at  fp. 

The  individual  spectral  points  used  in  computing  mean  spectra  are 
plotted  for  Lake  Worth  in  Appendix  B.  The  distribution  of  the  spectral 
points  about  the  mean  for  each  band  is  decidedly  unsymmetric  in  most 
cases.  The  distribution  usually  has  a  high  positive  skewness  indicating 
some  relatively  high  values  above  the  mean  in  the  distribution  which  are 
not  offset  by  correspondingly  low  values  below  the  mean.  However,  the 
distribution  of  spectral  points  about  the  mean  for  the  band  of  maximum 
energy  density  tends  to  be  more  symmetric  and  more  nearly  Gaussian. 

Third  and  fourth  moments  of  the  normalized  distribution  of  spectral 
points  in  selected  bands  were  computed  for  all  locations  considered  in 
this  report.  The  moments  are  defined  analogous  to  the  moments  of  the 
distribution  function  of  sea-surface  elevations  (eq.  6).  For  all  loca¬ 
tions,  both  the  third  and  fourth  moments  are  generally  lowest  in  the 
vicinity  of  the  spectral  peak  indicating  a  relatively  symmetric  and 
broad  distribution  about  the  mean. 

Typical  variations  of  the  third  and  fourth  moments  of  the  normalized 
distribution  of  spectral  points  as  a  function  of  frequency  are  shown  in 
Figures  14  and  15  for  mean  spectra  with  peak  periods  between  7.2  and  7.8 
seconds  at  Nags  Head.  The  third  and  fourth  moments  in  the  band  of  peak 
energy  density  actually  fall  below  the  values  for  a  normalized  Gaussian 
distribution.  Deviations  from  the  normalized  Gaussian  distribution  over 
all  frequency  bands  are  generally  less  for  high-energy  spectra  than  low- 
energy  spectra.  It  may  be  concluded  that  the  normalized  Gaussian  dis¬ 
tribution  gives  a  reasonable  representation  of  the  scatter  of  individual 
spectral  energy  densities  about  the  peak  of  each  mean  spectrum  in  Appen¬ 
dix  B,  but  it  is  not  reliable  for  representing  scatter  about  lower  energy 
bands  away  from  the  peak. 

Seasonal  variations  in  the  mean  spectrum  for  each  height-period  in¬ 
terval  were  studied  for  Nags  Head.  Four  3-month  seasons  were  selected 
using  the  seasonal  summaries  in  Thompson  (1977)  as  a  guide.  Mean  spectra 
were  confuted  for  winter  (December,  January,  and  February)  and  for  summer 
(June,  July,  and  August).  Seasonal  differences  in  the  mean  spectra  did 
not  appear  to  be  significant.  The  number  of  cases  per  height-period 
interval  was  considerably  reduced  when  the  spectra  were  grouped  by  sea¬ 
son.  Hence,  extensive  seasonal  summarization  was  not  considered  in  this 
study. 
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Figure  14.  Third  moment  of  the  normalized  distribution 
of  spectral  points  in  each  frequency  band 
as  a  function  of  frequency  and  H„  for  Nags 
Head,  North  Carolina  (T  «  7.2  to  7.8  seconds). 


Comparison  of  mean  spectra  in  Appendix  B  for  the  same  height-period 
intervals  at  different  locations  provides  additional  perspective.  Mean 
spectra  for  the  Atlantic  coast  locations  are  generally  quite  similar  to 
each  other  (Fig.  16, a  and  b).  Mean  spectra  for  the  gulf  coast  location 
differ  from  the  Atlantic  coast  spectra  in  that  spectral  energy  at  fre¬ 
quencies  lower  than  0.1  hertz  is  absent  in  virtually  all  cases.  Mean 
spectra  for  the  Great  Lakes  locations  show  an  even  more  pronounced 
absence  of  low-frequency  energy,  especially  for  low  wave  conditions. 

This  effect  is  illustrated  in  Figure  16a.  The  absence  of  low-frequency 
energy  gives  the  gulf  and  Great  Lakes  mean  spectra  a  tendency  to  have 
higher  peaks  than  the  mean  spectra  for  corresponding  height-period 
groups  at  Atlantic  coast  locations. 

Mean  spectra  for  the  Pacific  coast  locations  show  that  peak  periods 
longer  than  10  seconds  are  more  common  than  for  Atlantic  coast  locations. 
Reduced  amounts  of  energy  are  also  shown  at  frequencies  higher  than  about 
0.15  hertz  in  comparison  to  Atlantic  coast  spectra  with  long  peak  periods 
(Fig.  16b).  Hence,  the  Pacific  coast  long-period  mean  spectra  generally 
have  higher  peaks  than  their  counterparts  for  the  Atlantic  coast.  The 
few  Pacific  coast  mean  spectra  shown  with  peak  periods  shorter  than  9 
seconds  also  show  a  distinct  swell  energy  concentration  at  frequencies 
less  than  0.1  hertz. 

b.  Grouped  by  Height,  Period,  and  Water  Level.  Numerous  field 
studies  have  indicated  that  shallow-water  spectral  shape  is  dependent  on 
local  water  depth  (e.g. ,  Barnett  1969;  Hasselmann,  et  al.,  1973;  Wang  and 
Yang,  1976).  Therefore,  many  of  the  results  in  this  report  are  expected 
to  be  strongly  influenced  by  the  water  depth  at  the  gage  site.  The  higher 
wave  conditions  at  many  of  the  sites  can  represent  situations  where  waves 
occasionally  break  seaward  of  the  gage,  or,  in  the  extreme,  cases  where 
the  gage  is  actually  in  the  surf  zone  as  discussed  in  Section  II. 

To  give  more  perspective  on  the  systematic  influence  of  water  depth 
on  spectral  shape,  the  spectra  within  each  significant  height-peak  period 
interval  for  the  gage  at  Nags  Head  were  further  grouped  according  to  mean 
water  level  in  the  record.  Three  water  level  groups  were  arbitrarily 
defined.  Plots  of  average  spectra  by  water  level  for  height -period  in¬ 
tervals  in  which  there  were  a  reasonably  large  number  of  cases  are  shown 
in  Appendix  C.  Each  plot  is  labeled  with  the  height  and  period  intervals 
represented  and  with  a  mean  water  level  designator.  A  mean  equal  to  1 
indicates  the  mean  water  level  was  relatively  low;  a  mean  of  2  indicates 
midtide  water  levels  and  3  high  water  levels.  Water  level  2  encompassed 
a  range  of  49  centimeters.  The  differences  between  mean  spectra  as  a 
function  of  mean  water  level  are  small. 

3.  Spectral  and  Sea-Surface  Elevation  Distribution  Function  Parameters. 

a.  Secondary  Spectral  Peaks.  The  computerized  procedure  described  in 
Section  III  and  Appendix  D  was  used  to  identify  major  peaks  in  all  spec¬ 
tra  considered  for  this  report.  For  five  of  the  six  ocean  locations  only 
about  one-third  of  the  shallow-water  spectra  are  single  peaked  (Fig.  17). 
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Figure  17.  Percentage  of  spectra  versus 

number  of  major  spectral  peaks. 


At  the  remaining  ocean  location  (Pt.  Mugu)  only  5  percent  of  the  spectra 
were  single  peaked.  The  infrequency  of  single-peaked  spectra  at  Pt.  Mugu 
may  be  attributed  both  to  common  simultaneous  occurrence  of  sea  and  swell 
waves  and  to  pressure-compensation  effects  on  the  spectrum  which  often 
make  the  high-frequency  end  of  the  spectrum  look  like  a  meaningful  peak. 
At  all  ocean  locations,  about  one-half  the  spectra  are  double  peaked  with 
the  rest  containing  three,  four,  and  in  a  few  cases,  five  major  peaks. 
Several  spectra  with  five  major  peaks  and  corresponding  pen-and-ink 
traces  are  shown  in  Figure  18. 

Single-peaked  spectra  are  more  common  for  the  gulf  and  Great  Lakes 
sites  than  for  the  ocean  sites  (see  Fig.  17).  About  40  percent  of  the 
Naples  spectra  and  about  one-half  of  the  Great  Lakes  spectra  from  accel¬ 
erometer  buoys  6  kilometers  offshore  were  single  peaked.  Spectra  from 
the  Great  Lakes  nearshore  pressure  gages  show  that  only  20  to  30  percent 
of  the  cases  are  single  peaked.  However,  these  percentages  are  consid¬ 
ered  low  because  the  pressure-compensation  procedure  often  creates  a 
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Figure  18.  Pen-and-ink  strip-chart  records  and  corresponding  spectra 
for  some  cases  with  five  major  spectral  peaks:  (a)  Nags 
Head  record  at  0400  e.s.t.,  25  December  1968;  analysis 
at  0643  e.s.t.,  25  December  1968;  (b)  Nags  Head  record 
at  0800  e.s.t.,  24  February  1969;  analysis  at  0642  e.s.t., 
24  February  1969;  (c)  Huntington  Beach  record  at  1900 
e.s.t.,  9  May  1972;  analysis  at  1840  e.s.t.,  9  May  1972. 


small  peak  at  high  frequency.  Small  peaks  at  frequencies  less  than  0.05 
hertz  appeared  in  some  of  the  Great  Lakes  spectra,  especially  spectra 
from  the  Presque  Isle  pressure  gages  (App.  B) .  These  peaks  were  consid¬ 
ered  spurious  and  are  not  counted  as  major  spectral  peaks.  Double-peaked 
spectra  are  shown  to  be  common  at  both  the  gulf  and  Great  Lakes  locations. 


It  was  noted  that  secondary  spectral  peaks  do  not  always  represent 
independent  wave  trains.  When  waves  are  nonlinear,  i.e.,  when  H8/d  or 
He/Lp  are  large,  the  wave  profile  will  be  nonsinusoidal.  For  such  cases, 
a  secondary  spectral  peak  at  twice  the  dominant  frequency  is  expected 
even  though  a  photo  would  show  a  single  train  of  waves.  The  effect  is 
most  evident  for  high,  long  waves  at  low  tide;  i.e.,  cases  with  high 
values  of  Hg/d  and  low  values  of  d/Lp.  Appendixes  A  and  B  clearly 
show  that  secondary  peaks  tend  to  occur  systematically  for  high,  long 
shallow-water  waves.  However,  the  joint  distribution  table  of  number 
of  major  spectral  peaks  versus  significant  height  at  each  location  shows 
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no  clear  evidence  of  a  relationship  (see  Table  3).  The  correlation  be¬ 
tween  number  of  peaks  and  Hg  is  less  than  0.3  at  all  locations  (Table 
4).  The  joint  distribution  table  of  number  of  major  spectral  peaks  ver¬ 
sus  period  corresponding  to  the  highest  spectral  peak  at  each  location 
also  shows  little  evidence  of  a  relationship  (Table  5).  The  correlation 
between  number  of  peaks  and  Tp  is  also  very  low  (Table  4).  For  most 
wave  conditions  the  occurrence  of  secondary  spectral  peaks  does  not 
appear  to  be  systematically  related  to  Hg  and  Tp.  The  overall  sta¬ 
tistics  on  multiple  peaks  as  indicators  of  multiple  wave  trains  should 
be  reasonably  correct. 

Several  other  attempts  have  been  made  to  quantify  the  likelihood 
of  having  multiple  wave  trains.  Shepard  and  Inman  (1951)  report  visual 
wave  observations  to  identify  secondary  and  dominant  wave  trains  from 
the  Scripps  pier  in  southern  California  during  a  10-month  period.  They 
report  secondary  wave  trains  in  56  percent  of  302  observations.  This 
figure  may  be  an  underestimate  because  of  the  difficulty  of  visually 
identifying  secondary  trains,  but  is  reasonably  consistent  with  the 
Huntington  Beach  gage  data. 


Table  3.  Percentage  of  cases  in  each  significant  height  group 


corresponding  to  number  of  major  spectral  peaks. 
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Table  4.  Linear  correlation  between  number  of  major  spectral  peaks, 

significant  height  and  period  corresponding  to  the  highest  peak. 


Location 

Correlation 

Cases 

H e  vs. 

No.  of  peaks 

Tp  vs. 

No.  of  peaks 

Atlantic  City,  N.J. 

0.12 

0.07 

1,037 

Virginia  Beach,  Va. 

0.10 

0  10 

982 

Nags  Head,  N.C. 

0.16 

0.06 

945 

Lake  Worth,  Fla. 

0.02 

o.os 

982 

Naples,  Fla. 

0.17 

0.05 

866 

Michigan  City,  Ind.  (buoy) 

0.24 

0.17 

246 

Michigan  City,  Ind.  (pressure) 

0.26 

0.24 

127 

Presque  Isle,  Pa.  (buoy) 

0.09 

0. 18 

110 

Presque  Isle,  Pa.  (pressure) 

0.19 

0.32 

119 

Huntington  Beach,  Calif. 

0.01 

0.17 

746 

Pt.  Mugu,  Calif. 

0.02 

0. 16 

799 

Table  5.  Percentage  of  cases  in  each  peak  period  group 


corresponding  to  number  of  major  spectral  peaks. 
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17 

- 

6 
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21 

50 
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Hoffman  (1978)  reports  on  the  occurrence  of  multipeaked  spectra  from 
hindcasts  generated  for  a  grid  point  in  the  North  Atlantic  Ocean  by  the 
Spectral  Ocean  Wave  Model  of  the  U.S.  Navy  Fleet  Numerical  Weather  Cen¬ 
tral,  Monterey,  California.  The  criterion  for  accepting  a  secondary 
peak  was  that  its  peak  energy  be  at  least  30  percent  of  the  energy  in 
the  highest  peak  and  that  it  be  reasonably  separate  in  frequency  from 
the  main  peak.  About  25  percent  of  572  cases  considered  were  multi¬ 
peaked.  There  was  a  tendency  for  multipeaked  spectra  to  be  more  common 
during  low  wave  conditions  than  during  high  wave  conditions. 

It  is  of  some  engineering  interest  to  know  whether  a  spectrum  with 
a  certain  Tp  value  is  more  likely  to  contain  major  secondary  peaks 
than  a  spectrum  with  a  different  Tp  value.  Evidence  at  many  locations 
shows  that  spectra  associated  with  certain  ranges  of  Tp  do  indeed  ex¬ 
hibit  a  preference  for  multiple  peaks.  At  Atlantic  City,  Virginia  Beach, 
and  Nags  Head  on  the  Atlantic  coast  and  at  Huntington  Beach  on  the  Pacific 
coast,  multipeaked  spectra  are  more  common  for  long  and  short  Tp  cases 
than  for  intermediate  Tp  cases  (shown  in  Fig.  19).  Multipeaked  spectra 
are  least  common  for  Tp  values  between  about  8  and  10  seconds  at  the 
three  Atlantic  coast  locations.  At  Huntington  Beach  multipeaked  spectra 
are  least  common  for  Tp  values  of  12  to  13  seconds.  At  Lake  Worth  and 
Naples,  there  is  no  clear  evidence  of  a  relationship  between  Tp  and 
the  occurrence  of  secondary  spectral  peaks.  Data  on  secondary  peaks  from 
Pt.  Mugu  and  the  Great  Lakes  locations  were  not  considered  because  of 
difficulties  with  the  compensated  pressure  spectrum  and  because  of  the 
small  number  of  Great  Lakes  cases. 

Additional  insight  on  the  occurrence  of  secondary  spectral  peaks  can 
be  obtained  from  a  detailed  look  at  their  characteristics  for  a  repre¬ 
sentative  Atlantic  coast  location  and  the  single  Pacific  coast  location. 

For  example,  Figure  19  shows  that  70  percent  of  the  Nags  Head  spectra 
with  Tp  between  6  and  7  seconds  have  major  secondary  peaks;  Figure  20(a) 
shows  in  a  histogram  that  the  periods  of  those  secondary  peaks  range  from 
2  to  17  seconds.  The  figure  also  shows  that  secondary  peak  periods  are 
fairly  evenly  distributed  over  the  range  of  3  to  17  seconds.  The  rela¬ 
tively  low  percentage  of  secondary  peaks  at  7  to  8  seconds  shown  in  Fig¬ 
ure  20(a)  is  artificially  low  because  the  CERC  analysis  procedure  gives  a 
nonuniform  resolution  in  wave  period  which,  tends  to  decrease  the  number 
of  cases  in  the  7-  to  8-second  intervals  and  increase  the  number  of  cases 
in  the  6-  to  7-  and  8-  to  9-second  intervals.  The  procedure  also  precludes 
any  cases  with  intervals  of  11  to  12,  13  to  14,  15  to  16,  and  17  to  20 
seconds. 

Although  the  histograms  in  Figures  20  and  21  show  how  often  secondary 
peaks  occur,  no  indication  is  given  of  how  much  energy  is  contained  in 
the  secondary  peaks.  Therefore,  another  curve  was  added  to  the  figures 
to  show  the  average  energy  for  secondary  spectral  peaks  in  each  period 
interval.  The  energy  associated  with  each  peak  was  estimated  by  the  pro¬ 
cedures  discussed  in  Section  III  and  shown  in  Figure  7.  Relatively  high- 
energy  secondary  peaks  tend  to  occur  in  Figure  20(a)  at  periods  slightly 
less  than  Tp  and  at  periods  of  12  to  14  seconds.  The  long-period  sec¬ 
ondary  peaks  generally  represent  swell  waves. 
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Fraction  of  Spectro  Confoininq  at  t  .'last 
One  Major  Secondary  Peak 


Figure  19.  Fraction  of  spectra  with  at  least  one  major 
secondary  peak  versus  Tp. 
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Figure  20.  Frequency  of  occurrence  and  average  energy  for 
secondary  peaks  as  a  function  of  secondary  peak 
period  at  Nags  Head,  North  Carolina. 
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Figure  20(b)  shows  characteristics  of  secondary  spectral  peaks  at  Nags 
Head  when  period  of  the  dominant  peak  is  12  to  13  seconds.  Secondary 
peak  periods  of  4  to  9  seconds  are  common.  In  most  cases  these  secondary 
peaks  represent  locally  generated  wave  energy  or  swell  waves  generated 
near  the  Nags  Head  coast.  The  average  energy  associated  with  these  peaks 
is  relatively  high.  Secondary  peaks  at  periods  longer  than  10  seconds 
are  almost  totally  absent. 

At  Huntington  Beach,  spectra  with  relatively  short  Tp  have  a  very 
strong  tendency  for  a  secondary  peak  at  long  period  (see  Fig.  21, a). 
Locally  generated  waves  at  Huntington  Beach  are  usually  superimposed  on 
long-period  swell  waves.  Huntington  Beach  spectra  with  long  Tp  have 
a  small  tendency  for  secondary  peaks  with  short  periods  (see  Fig.  21,b). 
They  also  have  a  small  tendency  for  secondary  peaks  at  long  period,  which 
presumably  correspond  to  secondary  swell  waves  superimposed  on  the  domi¬ 
nant  swell.  However,  most  Huntington  Beach  spectra  with  long  Tp  appear 
to  be  dominated  by  the  main  spectral  peak. 

b.  Spectral  Peakedness.  The  spectral-peakedness  parameter,  Q^, 
defined  by  equation  (4)  was  computed  for  each  spectrum.  The  spectral- 
peakedness  parameter  is  shown  on  each  plot  in  Appendix  A. 

Cumulative  distribution  functions  of  Qp  are  shown  in  Figure  22 
for  the  selected  locations.  The  Qp  values  cover  a  wide  range  at  most 
locations— from  less  than  one  to  greater  than  five.  Qp  values  tend  to 
be  smallest  for  the  Atlantic  coast  sites.  The  highest  Qp  values  tend 
to  occur  for  the  Great  Lakes  gages,  with  the  pressure  gages  having  higher 
Qp  values  than  the  buoy  gages.  Qp  values  for  the  gulf  and  Pacific  coast 
gages  tend  to  be  intermediate  to  values  for  the  Great  Lakes  and  Atlantic 
coast  gages  except  for  the  highest  5  percent  of  the  Pt.  Mugu  Qp  values. 
The  cumulative  distribution  function  for  Pt.  Mugu  shows  a  rapid  increase 
in  Qp  values  associated  with  probabilities  of  less  than  5  percent. 

This  feature  is  a  result  of  occasional  large  overcompensation  of  the 
pressure  spectrum  which  creates  a  high  spectral  peak  at  high  frequency. 
Several  cases  in  which  Qp  was  greater  than  nine  are  omitted  from  Fig¬ 
ure  22.  These  cases  all  correspond  to  very  low  wave  conditions  except 
for  several  overcompensated  Pt.  Mugu  cases. 

Goda  (1976)  indicated  that  high  values  of  Qp  are  associated 
with  strong  grouping  of  high  waves.  Pen-and-ink  record  traces  for  sev¬ 
eral  cases  with  reasonably  high  waves  and  high  Qp  values  are  shown  in 
Figure  23.  Some  grouping  of  the  high  waves  is  evident  in  the  records. 
Records  with  low  Qp  values  in  Figures  13  and  18  show  little  evidence 
of  wave  grouping;  yet  the  record  in  Figure  11  for  high,  long  waves  shows 
evidence  of  high  waves  occurring  together,  but  the  Qp  value  is  only 
1.35. 

The  spectral-peakedness  parameter  might  be  expected  to  weakly  corre¬ 
late  with  the  number  of  major  spectral  peaks  so  that  single-peaked  spec¬ 
tra  tend  to  have  high  peakedness.  The  data  indicate  a  low  correlation 
at  some  locations,  but  ndt  at  others  (Table  6).  Qp  also  appears  to  be 
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Figure  22.  Cumulative  probability  distribution  functions 
for  spectral  peakedness. 
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Figure  23.  Pen-and-ink  strip-chart  records  for  selected  cases  with  high 
Qp  values:  (a)  Huntington  Beach  record  at  0100  e.s.t.,  3 
September  1972;  analysis  at  0040  e.s.t.,  3  September  1972; 

(b)  Naples  record  at  1220  e.s.t.,  9  June  1969;  analysis  at 
1340  e.s.t.,  9  June  1969;  (c)  Michigan  City  (buoy)  record 

at  2025  e.s.t.,  14  September  1975;  analysis  at  2220  e.s.t., 
14  September  1975. 
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Table  6.  Correlation  between  spectral-peakedness  parameter  and  other  wave  parameters. 


Location 

Cases 

Correlations 

Cases 

Correlations 

Tp  vs.  Qp 

vs. 

No.  of  peaks 

H 3  vs.  Qp 

(single 

-peaked  cases) 

Atlantic  City,  N.J. 

1,037 

0.22 

0.08 

322 

0.27 

Virginia  Beach,  Va. 

982 

0. 13 

0.02 

292 

0.20 

Nags  Head,  N.C. 

94S 

0.31 

0.20 

286 

0.20 

Lake  Worth,  Fla. 

982 

0.05 

0.08 

332 

0.29 

Naples,  Fla. 

866 

0.06 

0.16 

351 

0.31 

Michigan  City,  Ind.  (buoy) 

246 

0.09 

0.27 

128 

0.34 

Michigan  City,  Ind.  (pressure) 

127 

0.06 

0.61 

39 

0.80 

Presque  Isle,  Pa.  (buoy) 

110 

0.10 

0.  14 

54 

0.25 

Presque  Isle,  Pa.  (pressure) 

119 

0.06 

0.36 

25 

0.43 

Huntington  Beach,  Calif. 

746 

0.34 

0.24 

274 

0.05 

Pt.  llugu,  Calif. 

799 

0.22 

0.25 

34 

0.73 

weakly  correlated  with  relative  water  depth  as  evidenced  by  a  small  cor¬ 
relation  between  peak  period  and  Qp  for  single-peaked  spectra  (Table  6). 
This  may  also  be  a  result  of  the  definition  of  Qp.  A  low  correlation 
between  significant  height  and  Qp  is  also  shown  in  the  table. 

The  spectral-peakedness  parameter  appears  to  be  related  to  signif¬ 
icant  wave  height  and  peak  spectral  period  in  a  complicated  way  (see 
Table  7  for  Nags  Head  and  Huntington  Beach).  High  mean  values  of  Qp 
occur  for  significant  wave  heights  and  peak  periods  near  the  annual  mean 
for  most  locations  (Table  7).  Mean  Qp  values  decrease  fairly  system¬ 
atically  for  Tp  increasingly  higher  than  the  annual  mean  Tp  in  most 
cases  (Table  7) .  Qp  values  have  a  tendency  to  increase  with  Hs  for 
Tp  shorter  than  the  mean  and  decrease  with  Hs  for  He  and  Tp 
greater  than  the  annual  means.  High  values  of  Hs  are  often  associ¬ 
ated  with  relatively  high  mean  Qp  for  short  Tp,  but  low  mean  Qp 
for  long  Tp.  The  latter  is  mostly  a  result  of  the  nonsinusoidal  pro¬ 
file  for  long  waves  in  shallow  water. 

On  the  basis  of  evidence  in  this  report  and  Goda  (1976),  it  seems 
reasonable  to  hypothesize  that  high  peakedness  parameters  are  indicative 
of  strong  wave  grouping,  but  low  peakedness  parameters  in  shallow  water 
are  not  necessarily  indicative  of  weak  grouping,  especially  when  the 
relative  water  depth  is  small. 

c.  Skewness  and  Kurtosis  of  Sea-Surface  Elevation  Distribution 
Function"  The  distribution  function  for  sea-surface  elevation  also  pro¬ 
vides  useful  insight  on  shallow-water  spectral  statistics.  Mean  values 
of  skewness  and  kurtosis  (defined  in  eq.  6)  for  the  selected  locations 
are  given  in  Table  8.  It  is  evident  that,  except  at  the  Great  Lakes 
sites,  the  shallow-water  distribution  functions  tend  to  have  more  high 
than  low  extremes  (high  skewness)  and  that  they  tend  to  be  more  sharply 
focused  than  the  normalized  Gaussian  distribution  (high  kurtosis).  These 
findings  are  consistent  with  cnoidal  wave  profiles  shown  in  Figure  2. 


Cumulative  distribution  curves  for  skewness  and  kurtosis  are  shown 
in  Figures  24,  25,  and  26.  The  Atlantic  coast  locations  have  a  tendency 
for  higher  skewness  and  kurtosis  values  than  the  gulf  and  Pacific  coast 
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=  94  centimeters  and  Tr  *  8.9  seconds  for  Nags  Head; 


Table  8.  Statistics  on  skewr.ess  and  kurtosis  of  distribution  of  sea-surface  elevations 


Atlantic  City,  N.J. 

Virginia  Beach,  Va. 

Nags  Head,  N.C. 

Lake  North,  Fla. 

Naples,  Fla. 

Michigan  City,  Ind.  (buoy) 
Michigan  City,  Ind.  (pressur 
Presque  Isle,  Pa.  (buoy) 
Presque  Isle,  Pa.  (pressure) 
Huntington  Beach,  Calif. 

Pt.  Mugu,  Calif. 
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Figure  24.  Cumulative  probability  distribution  functions  for  skewness 
of  sea-surface  elevation  distribution  function. 
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Figure  25.  Cumulative  probability  distribution  functions  for  kurtosis 
of  sea-surface  elevation  distribution  function  at  Atlantic 
and  gulf  coast  locations. 
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Figure  26.  Cumulative  probability  distribution  functions  for  kurtosis 
of  sea-surface  elevation  distribution  function  at  Great 
Lakes  and  Pacific  coast  locations. 


locations.  The  Great  Lakes  locations  and  Pt.  Mugu  tend  to  have  the 
lowest  skewness  values.  The  Michigan  City  locations  and  the  Presque 
Isle  buoy  gage  show  a  tendency  for  occasional  very  high  kurtosis  values 
which  are  associated  with  very  low  significant  wave  heights.  Although 
several  kurtosis  values  greater  than  nine  occurred  for  the  Michigan  City 
locations  and  Pt.  Mugu,  the  values  are  omitted  from  Figures  25  and  26. 
The  Presque  Isle  pressure  gage  shows  a  tendency  for  exceptionally  low 
kurtosis  values.  A  section  of  pen-and-ink  record  for  a  case  with  high 
kurtosis  (equal  to  8.0)  was  shown  in  Figure  11. 

Some  of  the  kurtosis  values  in  Figures  25  and  26  are  questionably 
high.  Since  kurtosis  is  based  on  the  fourth  power  of  the  deviation  of 
sea-surface  elevations  from  the  mean  (i.e.,  n  =  4  in  eq.  6),  it  is  more 
affected  by  a  few  elevations  far  from  the  mean  than  any  other  parameters 
considered  in  this  report.  Noise  or  momentary  signal  losses  in  a  record 
which  are  comparable  to  the  size  of  the  largest  waves  will  pass  uncor¬ 
rected  through  the  editing  process  described  earlier.  These  bad  data 
points  usually  have  little  effect  on  the  spectrum  but  they  can  have  a 
significant  effect  on  the  kurtosis  of  the  distribution  of  sea-surface 
elevations. 

Another  factor  which  leads  to  questionably  high  values  of  kurtosis 
is  the  occasional  tendency  of  the  step-resistance  gage  to  become  elec¬ 
trically  shorted  at  a  particular  elevation,  often  by  marine  growth,  and 
to  produce  wave  records  in  which  the  troughs  are  unrealistically  flat. 
Because  kurtosis  is  especially  sensitive  to  these  problems  (known  to 
have  occurred  intermittently  in  the  data  considered  in  this  report). 


the  percentages  associated  with  kurtosis  values  in  Figures  25  and  26 
higher  than  about  4  should  be  considered  somewhat  inflated. 

The  correlation  between  skewness  and  peak  spectral  period  is  less 
than  0.25  at  all  locations;  for  single-peaked  cases  it  is  similarly  low. 
However,  the  correlation  between  skewness  and  significant  height  is 
moderatly  high  for  the  ocean  and  gulf  coast  locations  (Table  8).  There 
is  a  clear  tendency  for  skewness  to  increase  with  significant  height  for 
the  ocean  and  gulf  coast  locations  as  evidenced  in  a  plot  of  mean  skew¬ 
ness  versus  significant  height  (Fig.  27).  This  tendency  is  not  evident 
for  the  Great  Lakes  locations  (Fig.  28) . 

It  is  apparent  that  cnoidal  wave  profiles  (see  Fig.  2)  with  very 
high  skewness  and  kurtosis  values  are  rarely  observed  at  the  selected 
gage  locations.  However,  some  of  the  records  contain  wave  profiles  which 
are  more  nearly  cnoidal  than  sinusoidal.  Figures  24,  25,  and  26  indi¬ 
cate  that  cnoidal  profiles  appear  in  about  10  percent  of  the  records  or 
less. 

V.  CONCLUSIONS 

The  following  conclusions  are  drawn  from  evidence  and  discussion 
presented  in  this  report : 

1.  About  60  to  70  percent  of  the  ocean  and  gulf  coast  energy  spectra 
in  shallow  water  have  one  or  more  major  peaks  in  addition  to  the  dominant 
peak. 

2.  At  least  50  percent  of  the  spectra  in  Lakes  Michigan  and  Erie 
have  one  or  more  major  secondary  peaks. 

3.  Major  secondary  spectral  peaks  are  common  for  both  high  and  low 
energy  spectra. 

4.  The  very  small  sample  of  spectra  for  hurricane  waves  examined 
(seven  spectra)  does  not  reveal  any  special  characteristics  of  hurricane¬ 
generated  spectra. 

5.  Major  secondary  spectral  peaks  are  usually  indicative  of  secondary 
wave  trains. 

6.  When  secondary  spectral  peaks  occur  at  frequencies  of  two  or  three 
times  the  dominant  frequency,  they  may  represent  a  single  train  of  steep 
waves  rather  than  independent  wave  trains,  especially  if  the  relative 
water  depth  for  the  dominant  wave  is  small. 

7.  When  secondary  spectral  peaks  occur  at  frequencies  of  0.5  or  1.5 
times  the  dominant  frequency  and  the  waves  are  high  enough  for  breaking 
to  occur,  the  secondary  peaks  may  not  represent  independent  wave  trains. 


8.  There  is  evidence  that  wave  breaking  limits  the  ratio  Hs/d  at 
the  shallow-water  gage  locations  considered  to  less  than  0.55  rather 
than  the  commonly  used  limit  of  0.78. 

9.  Individual  shallow-water  ocean  spectra  with  comparable  signifi¬ 
cant  heights  and  peak  periods  are  highly  variable  in  shape. 

10.  Mean  spectra  for  the  same  significant  height-peak  period  group 
show  a  surprising  similarity  between  locations;  however,  there  is  a  nota¬ 
ble  lack  of  low-frequency  energy  in  gulf  and  Great  Lakes  spectra  and  a 
lack  of  high-frequency  energy  in  most  southern  Pacific  coast  spectra. 

11.  The  distribution  of  individual  energy  values  about  the  mean 
energy  for  each  spectral  band  tends  to  be  nearly  Gaussian  in  the  vicinity 
of  the  dominant  peak  but  non-Gaussian  for  low-energy  bands  away  from  the 
peak. 

12.  For  most  exposed  ocean  locations,  a  spectrum  is  more  likely  to 
contain  major  secondary  peaks  when  the  period  of  the  dominant  peak  is 
long  or  short  than  when  it  has  intermediate  values. 

13.  High  values  of  the  spectral-peakedness  parameter,  Qp,  indicate 
that  high  waves  may  be  occurring  in  groups;  however,  high  Qp  values  in 
shallow  water  are  usually  associated  with  low-energy  wave  conditions. 

14.  Low  values  of  Q^  generally  indicate  that  high  waves  are  poorly 

grouped.  However,  if  the  spectrum  has  secondary  peaks  which  do  not  repre¬ 
sent  independent  wave  trains,  is  not  a  direct  indicator  of  the  extent 

of  wave  grouping. 

15.  Qp  values  tend  to  be  higher  in  the  Great  Lakes  than  along  the 
ocean  and  gulf  coasts. 

16.  The  distribution  function  for  sea-surface  elevations  in  shallow 
water  tends  to  have  more  high  than  low  values  far  from  the  mean  and  tends 
to  show  less  variability  than  the  normalized  Gaussian  distribution. 

These  features  are  consistent  with  cnoidal-type,  shallow-water  wave 
profiles. 

VI.  SUMMARY 

Digital  wave  records  were  obtained  from  11  gages  at  9  shallow-water 
sites  along  the  U.S.  Atlantic,  Pacific,  gulf,  and  Great  Lakes  coasts. 

The  records  were  analyzed  with  a  fast  Fourier  transform  procedure  to 
produce  energy  spectra.  Selected  spectral  parameters,  including  the 
number  of  major  spectral  peaks  and  Goda's  (1970)  spectral-peakedness 
parameter,  and  parameters  of  the  distribution  function  of  sea-surface 
elevations  were  also  computed.  The  physical  meaning  of  the  shallow- 
water  spectrum  and  computed  parameters  was  discussed  and  illustrated 
with  three  cases  of  cnoidal  wave  profiles. 
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Spectra  for  the  24  highest  energy  cases  from  each  gage  were  summa¬ 
rized  individually.  Many  of  the  spectra  had  more  than  one  major  peak. 
Although  many  of  the  high-energy  cases  represent  breaking  wave  condi¬ 
tions,  the  ratio  of  significant  height  to  water  depth  was  less  than  0.55 
in  all  cases.  This  value  differs  significantly  from  the  value  0.78  which 
is  usually  used  in  conjunction  with  depth-limited  wave  heights. 

All  spectra  were  grouped  according  to  significant  height  and  period 
corresponding  to  the  highest  spectral  peak.  A  mean  spectrum  was  com¬ 
puted  for  most  height -period  groups.  The  standard  deviation  of  spectral 
values  about  the  mean  was  also  included.  The  mean  spectra  show  a  clear 
lack  of  low-frequency  energy  at  the  gulf  and  Great  Lakes  sites  and  a 
lack  of  high-frequency  energy  at  the  southern  California  site.  They  also 
show  evidence  of  systematic  variations  in  shape  as  a  function  of  signifi¬ 
cant  height  and  peak  period. 

Multipeaked  spectra  are  common;  60  to  70  percent  of  the  ocean  and 
gulf  coast  c  ectra  and  at  least  50  percent  of  the  Great  Lakes  spectra 
have  more  cnan  one  major  peak.  Major  spectral  peaks  are  usually  indica¬ 
tive  of  independent  wave  trains. 

High  values  of  the  spectral-peakedness  parameter  appear  to  be  related 
to  the  occurrence  of  high  waves  in  groups.  However,  low  values  of  the 
spectral-peakedness  parameter  are  not  necessarily  indicative  of  a  lack 
of  grouping  of  high  waves  in  shallow  water.  The  distribution  function  of 
sea-surface  elevations  in  shallow  water  tends  to  deviate  from  a  normalized 
Gaussian  distribution  in  a  way  which  is  consistent  with  cnoidal-type  wave 
profiles . 

In  the  course  of  this  investigation,  a  computer  program  was  developed 
to  identify  major  peaks  and  valleys  in  an  irregular  signal.  The  program 
has  many  potential  applications,  including  identification  of  major  peaks 
in  a  spectrum  and  identification  of  meaningful  crests  and  troughs  in  a 
time  series  of  sea-surface  elevation. 
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INDIVIDUAL  HIGH-ENERGY  WAVE  SPECTRAL  PLOTS 
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energy  wave  spectral  plots,  Atlant 


gure  A- 1.  High-energy  wave  spectral  plots,  Atlantic  City,  New  Jersey — Continued 
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energy  wave  spectral  plots,  Virginia  Beach,  Virginia. 


energy  wave  spectral  plots,  Virginia  Beach,  Virginia. — Continued 
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energy  wave  spectral  plots,  Nags  Head,  North  Carolina 
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Figure 


energy  wave  spectral  plots.  Lake  Worth,  Florida. --Continued 
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energy  wave  spectral  plots,  Naples,  Florida. 


energy  wave  spectral  plots,  Naples,  Florida. --Continued 


17  Feb.  0141  91  4.52  S.47  0.2012  0.167  0.0336  0.0273  1.63  0.487  3.48 
17  Aug.  1541  89  8.79  5.89  0.0929  0.152  0.0141  0.0078  1.35  0.407  3.33 
9  Feb.  0941  83  5.97  5.23  0.1358  0.169  0.0229  0.0150  2.20  0.468  3.55 
4  Mar.  1542  88  4.31  5.19  0.2074  0.169  0.0350  0.0285  1.77  0.322  3.48 


Figure  A-6.  High-energy  wave  spectral  plots,  Michigan  City,  Indiana  (buoy). 
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re  A-6.  High-energy  wave  spectral  plots,  Michigan  City,  Indiana  (buoy) . --Continued 
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Figure  A-7.  High-energy  wave  spectral  plots,  Michigan  City,  Indiana  (pressure  gage) . --Continued 
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Figure  A-10.  High-energy  wave  spectral  plots,  Huntington  Beach,  California. — Continued 
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Table  B-l.  Approximate  relative  depth  and  ratio  of  signif¬ 
icant  height  to  depth  for  mean  spectra  grouped 
by  height  and  period. 


|  Parameters  based  on  period 

|  Parameters  based  on  height 

KB 

d/Lp 

»6 

He/d 

(cm) 

r 

Atlantic  City  | 

wmwm 

iMnifE  3B 

46 

■m 

witrm 

76 

ns 

0.149 

107 

HIS 

0.124 

0.013 

137 

6.8 

0.115 

0.011 

168 

0.324 

7.4 

0.105 

0.010 

198 

0.382 

■m 

0.096 

0.008 

229 

0.441 

8.8 

0.086 

0.007 

259 

0.500 

9.7 

0.078 

0.006 

10.5 

0.069 

0.005 

12.3 

0.060 

0.003 

14.1 

0.052 

0.003 

16.7 

0.044 

0.002 

|  Virginia  Beach  | 

3.5 

0.327 

0.051 

46 

0.075 

4.5 

0.218 

0.031 

76 

0.125 

5.5 

0.165 

0.021 

107 

0.175 

6.4 

0.136 

0.015 

168 

0.275 

6.8 

0.127 

0.013 

198 

0.325 

7.4 

0.115 

0.011 

229 

0.375 

8.0 

0.105 

0.010 

259 

0.425 

8.8 

0.094 

0.008 

290 

0.475 

9.7 

0.085 

0.007 

10.8 

0.075 

0.005 

12.3 

0.066 

0.004 

14.1 

0.057 

0.003 

|  Nags  Head  | 

5.S 

0.149 

0.017 

46 

0.088 

6.4 

0.124 

0.013 

76 

0.147 

6.8 

0.115 

0.011 

107 

0.206 

7.4 

0.105 

0.010 

137 

0.265 

8.0 

0.096 

0.008 

168 

0.324 

8.8 

0.086 

0.007 

198 

0.382 

9.7 

0.078 

0.006 

229 

0.441 

10.8 

0.069 

0.005 

259 

0.500 

12.3 

0.060 

0.003 

14.1 

0.052 

0.003 

16.7 

0.044 

0.002 
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Table  B-l.  Approximate  relative  depth  and  ratio  of  signif¬ 
icant  height  to  depth  for  mean  spectra  grouped 
by  height  and  period. — Continued  _ 


|  Parameters  based  on  period 

|  Parameters  based  on  height  ] 

TP 

d/1 v 

d/gT| 

Hs 

Hg/d 

(s) 

(cm) 

|  Lake  WoTth  | 

2.5 

0.563 

0.089 

46 

0.083 

3.5 

0.299 

0.046 

76 

0.139 

4.5 

0.202 

0.028 

107 

0.194 

5.5  , 

0.154 

0.018 

137 

0.250 

6.4 

0.128 

0.014 

168 

0.306 

8.8 

0.007 

198 

0.361 

9.7 

10.8 

KIR  tfl 

12.3 

14.1 

imjtzM 

|  Naples  ] 

■a 

15 

■sis 

46 

hi 

WSSm 

76 

0.149 

WSSmv 

107 

0.124 

137 

6.8 

0.115 

0.011 

7.4 

0.105 

0.010 

8.0 

0.096 

0.008 

8.8 

0.086 

0.007 

9.7 

0.006 

Michigan  City  (buoy) 


1 

| 

0.298 

0.152 

0.092 

0.062 

0.045 

0.040 

0.034 

0.029 

15 

46 

76 

107 

137 

168 

198 

229 

259 

290 

0.008 

0.025 

0.042 

0.058 

0.075 

0.092 

0.108 

0.125 

0.142 

0.158 

Michigan  City  (pressure)  | 

3.5 

15 

4.5 

B  I 

0.035 

46 

0.065 

5.5 

■m 

0.024 

76 

0.109 

6.4 

0.149 

0.017 

107 

0.152 

6.8 

0.138 

0.015 

137 

0.196 

7.4 

0. 12S 

0.013 

168 

0.239  ( 

8.0 

0.114 

0.011 

198 

0.283 

229 

0.326 

259 

0.370 
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Table  B-l.  Approximate  relative  depth  and  ratio  of  signif¬ 
icant  height  to  depth  for  mean  spectra  grouped 
_ by  height  and  period. —Continued 


Paramct 

ers  based  on  period 

j  Parameters  based  on  height 

Tf> 

d/1 v 

d/gTp 

!  H* 

Ha/d 

(s) 

1  («") 

Presque  Isle  (buoy) 

3.5 

0.973 

0.155 

46 

0.025 

4.S 

0.588 

0.094 

76 

0.041 

5.5 

0.397 

0.063 

107 

0.057 

6.4 

0.301 

0.046 

137 

0.074 

6.8 

0.272 

0.041 

168 

0.090 

7.4 

0.238 

0.035 

198 

0.107 

229 

0.124 

Presque  Isle  (pressure)  1 

3.5 

0.257 

0.038 

15 

0.033 

4.5 

0.178 

0.023 

46 

0.100 

!  5.5 

0.138 

0.015 

76 

0.167 

6.4 

0.115 

0.011 

107 

0.233 

6.8 

0.107 

0.010 

137 

0.300 

7.4 

0.097 

0.009 

168 

0.367 

1  Huntington  Beach 

5. 5 

0.219 

0.031 

46 

0515 

6.4 

0.023 

76 

0.083 

6.8 

0.020 

107 

7.4 

0.017 

138 

cTSVHSEi 

8.0 

0.015 

168 

8.8 

■  I 

0.012 

198 

9.7 

0. 106 

0.010 

229 

10.8 

0.094 

0.008 

259 

0.283 

12.3 

0.006 

290 

0.317 

14.1 

0.005 

16.7 

0.059 

0.003 

20.3 

0.048 

0.002 

—  - 

Pt. 

Mugu 

3.5 

46 

0.048 

4.5 

■  ■ 

1 

76 

0.081 

5.5 

IKZIZXiil 

107 

0.113 

6.4 

0.181 

137 

0.145 

6.8 

0.166 

0.021 

168 

0.177 

7.4 

0.150 

0.018 

198 

0.210 

8.0 

0.136 

0.015 

229 

0.242 

8.8 

0.121 

0.012 

9.7 

0. 108 

0.010 

10.8 

0.096 

0.008 

12.3 

0.083 

0.006 

14.1 

0.072 

0.005 

16.7 

0.060 

0.003 

20.3 

0.049 

0.002 
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APPENDIX  C 

MEAN  SPECTRA  GROUPED  BY  HEIGHT,  PERIOD,  AND  WATER  LEVEL 


Interpretation  key 


Mean 

Water  level 
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low 

2 
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APPENDIX  D 

DESCRIPTION  OF  COMPUTER  ROUTINE  SMOOTH  FOR  IDENTIFYING 
MAJOR  PEAKS  AND  VALLEYS  IN  AN  IRREGULAR  SIGNAL 

The  computer  routine  SMOOTH  is  useful  for  deleting  small,  inconse¬ 
quential  peaks  and  valleys  from  an  irregular  digital  signal.  Peaks  and 
valleys  that  remain  after  application  of  SMOOTH  represent  major  extrema 
which  in  many  cases  are  more  meaningful  than  small  wiggles  in  the  signal. 

SMOOTH  is  a  very  general  routine  which  can  be  applied  to  many  differ¬ 
ent  situations.  It  has  been  applied  at  CERC  to  the  following  three 
situations: 

(a)  Confutation  of  statistics  on  individual  waves.  SMOOTH  was 
applied  to  a  time  series  of  sea-surface  elevations; 

(b)  Computation  of  statistics  on  tidal  highs  and  lows.  SMOOTH 
was  applied  to  a  time  series  of  tidal  elevations. 

(c)  Identification  of  major  spectral  peaks.  SMOOTH  was  applied 
to  a  wave  energy  spectrum. 

The  operation  of  SMOOTH  is  most  conveniently  described  in  terms  of 
its  application  to  a  time  series  of  sea-surface  elevations,  although  its 
other  applications  are  analogous.  The  general  scheme  of  operation  con¬ 
sists  of  a  check  on  the  time  difference  and  elevation  difference  between 
successive  extrema.  If  either  is  less  than  the  specified  acceptable  min¬ 
imum  then  one  peak  and  one  valley  are  deleted  from  the  time  series. 

The  input  to  SMOOTH  consists  of  several  control  parameters  and  an 
array  (EXTIM)  containing  time  and  elevation  for  each  extremum  in  the 
time  series.  Figure  D-l  shows  five  extrema  in  a  hypothetical  time 
series.  If  the  point  labeled  "-1"  were  the  first  point  in  the  time 
series,  then  the  first  10  values  in  the  EXTIM  array  would  be 

t_j,  0_i»  tQ,  Hq ,  tj,  rij,  t2»  D2»  t3,  and  113, 

where  and  n-£  are  defined  as  the  time  and  elevation  associated 

with  the  ith  point.  The  control  parameters  which  must  be  specified 
are: 

FURST  *  time  associated  with  the  first  point  desired  in  the 
time  series. 

ITEMS  =  total  number  of  values  in  EXTIM  array  (=  twice  the 
number  of  extrema) . 

CHP  *  minimum  acceptable  time  difference  between  successive 
points  (critical  half  period). 

HMIN  *  minimum  acceptable  elevation  difference  between 
successive  points. 


145 


'V31 


<vv 

Figure  D-l.  Hypothetical  time  series. 

After  completion  of  SMOOTH,  the  times  and  elevations  of  major  peaks 
and  valleys  are  stored  in  the  first  ITEMS  elements  of  EXTIM.  The  value 
of  ITEMS  has  been  reduced  in  accordance  with  the  number  of  small  peaks 
and  valleys  eliminated.  Since  the  smoothing  algorithm  cannot  work  prop¬ 
erly  at  the  end  of  a  record,  the  last  few  points  are  usually  accepted 
regardless  of  whether  or  not  they  satisfy  the  acceptance  criteria.  Thus, 
the  last  six  elements  in  EXTIM  may  not  represent  major  peaks  or  valleys 
and  should  either  be  checked  or  categorically  eliminated. 

A  complete  list  of  subroutine  SMOOTH  with  comments  is  provided  in 
Figure  D-2.  A  flow  chart  for  the  main  processing  loop  (identified  as 
loop  1520  in  the  listing)  is  provided  in  Figure  D-3.  The  symbols  n£ 
and  t^  and  the  point  numbers  referred  to  in  the  flow  chart  are  con¬ 
sistent  with  definitions  in  Figure  D-l.  ^CHP  tiMIN  in  the  ^low 

chart  are  equivalent  to  CHP  and  HMIN  in  the  listing. 
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SUBROUTINE  SMOOTH  (FU«$T,  ITEMS) 

C  SUBROUTINE  SHDUTm  T  *NF S  A  HE  CORO  OF  PFAKS  AND  WALLERS  AND 

C  ELIMINATES  TNCONS£UU{ NTIaL  PfakS  AND  VaLLFTS,  The  ELIMINATION 

C  CRITERIA  are  *  minimum  ojfffrfnCE  BETwEEn  PEAK  and  valley 

C  ELEVATIONS  (MMJNJ  AMO  A  minimum  horizontal  SPACING  (BE  IT 

C  TIME,  POSITION.  ETC)  PETmEEn  PEAK  AND  V*LLtV  CCMP), 

c  input  parameters  are  defined  as  follows, 

C  Fijas'T  c  STARTING  time 

C  ITEMS  ■  TOTAL  number  OF  EXTREME  VaLOCS  (INCLUDE*  TIME  AMO 
C  ELEVATION  VALUES) 

C  CMP  ■  CRITICAL  malt  PERIOD 

C  mp I n  ■  minimum  HEIGHT  TO  bE  CONStOEREO 

C  EXTIM(ODO)  *  TIME 

C  EXTIM(EVEN)  «  EXTREME  ELEVATION  ASSOCIATED  *XTH  EXTIM(evEN.j)  VALUE 

C  OF  TIRE 

common  /SMOT/  EXTJMTlOO) 

C  SET  VALUES  OF  CMP  AND  mmin  TO  BE  USEO 

DATA  CMP,  MMIN  /  O.OOi  J.O  / 

c  initialize  variables  and  find  starting  point  for  processing 

1503  IIT£Ms»IT£ms 
J«1 

ITEmT«1 TEmj«e 

DO  150?  !■) ,  I  ITEMS, 2 

ISAiP»I 

IFIEXTIMtn.GE.FURST)  GO  TO  1504 

JmJ*l 

I  TEMS* I TEMSmJ 
150?  CONTINUE 

1504  ISTART»ISXIP 

CMAMFMMMMMMMMTMMflUMMIMMMFMMMMMmtFA.MMFTMMM 

C  BEGIN  MAIN  PROCESSING  LOOP 

00  1520  I»ISTaRT.ITEM?,2 
L»I 

IF (ISTART.GT.ITem7)G0  TO  1520 
IF(tXTlM(l)tCMP,GT.t*TlHf I«?))GO  TO  151  I 
IF (AbSttXTIMf j4j )*EXTlM(lt3)),LT,HMlN)G0  TO  1511 
C  IF  no  TRANSFER,  TMIS  EXTREME  ACCEPTED 
Exti*, J)*EXTIM(J) 

ExTlM(J*i )bExTIm(1*1) 

JeJ*2 
GO  TO  1520 

C  RM£N  The  NEXT  instruction  is  RE^Cmed,  one  high  and  ONE  LO"  WILL 
C  BE  DELETED 

1511  IF(E«T!M(IM),GT.tXTI«<lAj))GO  TO  1514 

C  IF  NEXT  INSTRUCTION  1»  used,  TmIS  is  *  LOR 
IFfE*TIMClAB).GT.ExTIMUM))GO  TO  1512 
IFtEXTlM(I*l),GT,E*TIM{l»i))GO  TO  1515 
GU  TO  1517 

1512  IF(EXTIM(Nj),GT,ExTIM(it7jiGo  TO  1519 
GO  TO  151b 

C  when  NEXT  INSTRUCTION  Is  USEO,  THIS  IS  A  HJGH 
15l«  IF(EXTIM(I*5),GT,EXTIH(I*1))G0  TO  1515 
IF(EXTlMU*7),GT,fcXTlM<l»J))G0  TO  1519 
GO  TO  1516 

1 5 1 5  IF  (EXTIM(J»J),GT.EXTIM{I.1))GD  TO  »Sl7 
GO  TO  1518 

c  set  the  value  of  zcasc 

1516  ICASEal 

60  TO  1521 

1517  ICASER2 
GO  TO  1521 

Figure  D-2.  List  of  subroutine  SMOOTH. 
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1518  1CASE«3 
GO  TO  1521 
15l«  IC A3E*R 

1521  Jlal»2 
J2*I*» 

c  delete  one  high  anq  une  lo*  according  to  the  value  of  icase 

GO  TO(1522,1523,l52a,l52S)lCASE. 

C  storage  plan  a 

1522  £xTlH{J-2)atxTIh(I«2) 

EXTjH(J-n*EXTlHtl»l> 

E  X  T  JH( j jafcxt  JH( 1 ) 
t»U«lJ»l  )sEXTIh(1ai) 

ExTI«(J»2)stXT!H(lA6) 

EXTIMf J*l)*EXTlMfl»7) 

GO  TO  1526 

C  STORAGE  PLAN  B 

1521  EXTlM(J-2)aExTjN(I*2) 

E*TlH(J.j)aEXTlMtl.i) 

E*Tl-(  J)sfcxTIHU*«) 

EXTlH(JM)*tXTlHtI»5> 

ExTIH<J»2)afc*TlH(I*6) 

ExTI*( Ja3)*ExT1h(1*7) 

GO  TO  1526 

C  STORAGE  PLAN  C 

1526  E*TlHCJ-2)aEXTlM(I*2) 

ExT]"(J“1)*£xtIn(1aJ) 

ExT1h(j)sextim( !♦«) 

EXTJN{J*l)stXTlH(I*5) 

EXTI~( J*2)SE xTlHf ]66) 

ExTjm(  J*3)«exTIH{  1-675 
GO  TO  1526 
:  STORAGE  PLAN  0 

1525  ExTIh(I+6)sExTIH(I+2) 

ExTlN(lt7)8EXTlM(l63) 

ExTlH(J-2),tXTlHf 1*2) 

ExTIH( J.i)*txTlH(J*l) 

ExTlH{J)a£xTIH(l) 

ExTIH(Jti)atXTI“n*U 

ExTlH( jA2JatXT|H(lA2] 

EXTlH(JfjjatXTlH(lAJ) 

Kial*q 
52*  I  ♦  1  0 

1526  J 1  a J«2 
J2aJ*3 

ITE «S»IT£H$-c 
IS*IP*It8 

GO  TO  150fl 
1520  CONTINUE 

:  '  end  ha  In  PROCESSING  LOOP 

. . 

JxJ-1 

IMI»2.tO.!TEHT)L*l 
1530  IE (L.GT,1ITE“S)G0  TO  15«0 
Ja  J*  | 

ExTlHJ J)«tXTlH(L) 


L«Lt1 

GO  TO  1530 
lSuO  CONTINUE 

"ERE  ANT  OELFT1UNJ  IN  THIS  PASS  RtPROCESS  ALL  RtHA|NJNG 
hIg»S  ANO  cn.S  TO  HAKF  FURTHER  DELITIOnS  If  NEEDtD( 
IFdlTEHS.GT.ITEHSlGO  TO  150 J 
RETURN 

end 

Figure  D-2.  List  of  subroutine  SMOOTH  .--Continued 


Figure  D-3.  Flow  chart  for  main  processing  loop  in  SMOOTH  (loop  1520). 
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